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The solar neutrino capture rate measured by the Russian- American Gallium Experiment (SAGE) 
on metallic gallium during the period January 1990 through December 1997 is 67.2^y Q^3'Q SNU, 
where the uncertainties are statistical and systematic, respectively. This represents only about half 
of the predicted standard solar model rate of 129 SNU. All the experimental procedures, including 
extraction of germanium from gallium, counting of ^^Ge, and data analysis, are discussed in detail. 



PACS number(s): 26.65.+t, 95.85.Ry 13.15.+g 
I. INTRODUCTION 

The Sun produces its energy by the nuclear fusion of 
four protons into an a particle, chains of reactions that 
yield two positrons and two neutrinos. Since these low- 
energy neutrinos are weakly interacting, it was assumed 
that they traverse the Sun and reach the Earth without 
change. Measurement of the neutrino energy spectrum 
should thus give information about the conditions un- 
der which the nuclear reactions take place in the Sun. 
All solar neutrino experiments, however, have observed 
considerably fewer neutrinos than are predicted by de- 
tailed models of the physical processes in the Sun that 
are based on the nuclear reaction chains. As a result of 
this neutrino deficit, the assumption that the neutrinos 
are unchanged during their passage from the Sun to the 
Earth is now seriously questioned. For such transforma- 
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tions to occur neutrinos must have mass, a hypothesis of 
far-reaching consequences. 

The experimental study of solar neutrinos is now over 
30 years old. The first experiment, a radiochemical de- 
tector based on chlorine ||l ,^ , observed a capture rate of 
2.55±0.17±0.18 SNU, where 1 SNU = 1 interaction/s in 
a target that contains 10^^ atoms of the neutrino absorb- 
ing isotope. Although different standard solar models 
(SSM's) predict somewhat different rates for the chlorine 
experiment (for example, SNU |,| and 7.2 SNU 

Q), all such models predict a rate significantly higher 
than observed. 

For 20 years, until about 1985, the chlorine experiment 
was the only measurement. This experiment is primar- 
ily sensitive to high-energy ^B neutrinos with a ^ 20% 
contribution from other sources, mainly ^Be. The flux 
of ^B neutrinos is very dependent on the central tem- 
perature of the Sun {T^ As a result many mod- 
els were suggested that would slightly suppress and 
hence decrease the ®B flux significantly. (See Ref. M for 
a description of a large number of such models.) Most 
of these models, however, run into difficulty with some 
other measured aspect of the Sun. An alternative solu- 
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TABLE I. Predicted solar neutrino fluxes M and their contribution to the Ga capture rate. 





V 


V energy 


V flux 


Ga capture 


Reaction 


branch 


(MeV) 


(cm^^ s"-^) 


rate (SNU) 


p + p^d + 6+ + 


pp 


0-0.42 


(5.94 ±0.06) X 10'" 


69.6 


'^Be + e-^^B + iy 


^Be 


0.38, 0.86 


(4.80 ± 0.43) X 10'-^ 


34.4 


«B^«Be* + e+ + z/ 


«B 


0-14.1 


(5.15ir72) X 10' 


12.4 


CNO reactions 


CNO 


0-1.73 


(1.1 ±0.2) X 10** 


9.8 


p + e~ + p—>d + V 


pep 


1.44 


(1.39 ±0.01) X 10** 


2.8 



tion to this discrepancy could be neutrino oscillations. 
The CI experiment operates on the inverse beta decay 
reaction and thus is only sensitive to electron neutrinos. 
If the neutrinos were to change flavor on their trip from 
the solar core to the Earth, the CI experiment would not 
observe them. 

In the mid 1980s, the Kamioka nucleon decay experi- 
ment (Kamiokande) began to measure the solar neutrino 
flux. This large water Cherenkov detector was originally 
designed to look for high-energy signals from proton de- 
cay. After great effort, the energy threshold was reduced 
to a level to permit a sensitivity to recoil electrons from 
solar neutrino interactions. The path of the recoil 
electrons is in the direction of the initial neutrino tra- 
jectory, and thus this experiment demonstrated for the 
first time that neutrinos were coming from the Sun. The 
measured flux § of (2.80 ± 0.19 ± 0.33) x 10V(cm2 s) 
was less than half of the solar model prediction, and the 
solar neutrino problem was thus confirmed by a second 
experiment. 

Because the high-energy solar neutrino flux was sup- 
pressed, it became very important to also determine the 
flux of low-energy neutrinos produced in the dominant 
proton-proton (pp) reaction. Exotic hypotheses aside, 
the rate of the reaction is directly related to the solar 
luminosity and is insensitive to alterations in the solar 
model. In the early 1990s the Russian- American Gal- 
lium Experiment (SAGE) and then the Gallium Exper- 
iment (GALLEX) began to publish results. These ex- 
periments are based on the neutrino capture reaction 
7iGa(i/e,e-)^iGe § and have the very low threshold 
of 233 keV They are thus sensitive to low-energy 

pp neutrinos, whose end point energy is 423 keV 
and provide the only feasible means at present to mea- 
sure low-energy solar neutrinos. The SAGE result |lj] 
of 66.9+^ jt^;^ SNU with a target of Ga metal and the 
GALLEX result (l| of 77.5 ± 6.2^^;? SNU with a target 
of GaCla are both well below the SSM prediction from 
the Bahcall-Pinsonneault solar model Q of 129l6 ^NU. 
The insensitivity of Ga to the solar model is seen in the 
capture rate calculation from the model of Brun, Turck- 
Chieze, and Morel || of 127.2 SNU. The contributions 
of the components of the solar neutrino flux to the ''^Ga 
capture rate are given in Table ^ 

With the four experiments having three different 
thresholds, one can deduce some information concerning 
the V energy spectral distribution. If one fits the data 



from all experiments with the neutrino fluxes as free pa- 
rameters, the best fit is when the t'Bc flux is greatly re- 
duced whereas there is an appreciable flux [O-M. 
This is an apparent paradox as it is difhcult to form^3 
in the Sun without forming ''Be. 

In 1996 Super-Kamiokande began to take data. This 
50 kton water Cherenkov detector is the first high-count- 
rate solar neutrino experiment. The present result IT^ ] 
for the ^B neutrino flux, assuming that neutrino trans- 
formations do not occur, is (2.44±0.05iS5:E5?) x lOVlcm^ 
s), in agreement with its predecessor. 

The purpose of this paper is to summarize all of the 
SAGE data for the last eight years. It is organized in 
the same way as the SAGE experiment is carried out: 
after presenting some general aspects of the experiment, 
we consider the chemical extraction of Ge from Ga and 
the subsequent Ge purification. Then we present how the 
Ge is counted, how ^^Ge events are identified, and how 
the data are analyzed to give the solar neutrino capture 
rate. Finally, we consider the sources of systematic un- 
certainty, give the overall results, and conclude with the 
implications for solar and neutrino physics. 

In an attempt to make the material understandable to 
the general reader, but still useful to the specialist, each 
of these subjects is first discussed in a general way, fol- 
lowed by subsections that give more detail. The reader 
who wants a general overview need only read the begin- 
ning of each section. The reader who desires more in- 
formation regarding a particular subject should read the 
appropriate subsection. 

II. SAGE OVERVIEW 

In this section we give some general information on the 
location of the experiment, its physical characteristics, 
and the division of the SAGE data into three experimen- 
tal periods. 

A. Baksan Neutrino Observatory 

The SAGE experiment is situated in a specially built 
underground laboratory at the Baksan Neutrino Obser- 
vatory (BNO) of the Institute for Nuclear Research of the 
Russian Academy of Sciences in the northern Caucasus 
mountains. The main chamber of the laboratory is 60 m 
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TABLE II. 
data. 



Definition of the various segments of SAGE 



Designation Included extractions 



Comments 



SAGE I 
SAGE II.l 

SAGE II.2 
SAGE II.3 
SAGE III.l 



Jan. 90 
Sep. 92 

Nov. 93 
July 94 - 
Jan. 95 



May 92 
Oct. 93 

■ June 94 
Dec. 94 
June 95 



SAGE III.2 July 95 ^ present 



Rise time from ADP 
Rise time from wave form, 
begin to use L peak 

Ga theft period 
Before Cr experiment 

Some extractions 
during Cr experiment 
After Cr experiment 



long, 10 m wide, and 12 m high. It is located 3.5 km 
from the entrance of a horizontal adit excavated into the 
side of Mount Andyrchi and has an overhead shielding 
of 4700 meters of water equivalent. To reduce neutron 
and gamma backgrounds from the rock, the laboratory 
is entirely lined with 60 cm of low-radioactivity concrete 
with an outer 6 mm steel shell. All aspects of the experi- 
ment are in this underground area, with additional rooms 
devoted to chemistry, counting, and a low-background 
solid-state Ge detector. Other facilities for subsidiary 
measurements are in a general laboratory building out- 
side the adit. 



B. Extraction history 

The data from SAGE span nearly a decade during 
which the experiment evolved a great deal. As a result, 
the data can be naturally divided into several periods 
characterized by different experimental conditions. Ex- 
tractions on approximately 30 tons of Ga began in 1988; 
by late 1989 backgrounds were low enough to begin so- 
lar neutrino measurements. The data period referred to 
as SAGE I began in January 1990 and ended in May 
1992 In the summer of 1991, the extraction mass 

was increased to nearly 60 tons. The SAGE I data were 
taken without digitized wave forms and the L peak could 
not be analyzed because of high electronic noise at low 
energy. (The decay modes of ''-^Ge are described below 



in Sec. IV.) The solar neutrino capture rate determined 
from this data were published in Ref. ||l^ . 

Within a few months after SAGE I was completed, the 
experiment was greatly improved with respect to elec- 
tronic noise. The following period of data, from Septem- 
ber 1992 to December 1994, is referred to as SAGE II. It 
is distinguished by recording of the counter wave form in 
most runs which makes possible analysis of events in the 
low-energy L peak. 

During SAGE II, there was a period (which we call 
SAGE II. 2) in which 2 tons of gallium, approximately 
3.6% of the total mass, was stolen from the experiment. 
The gallium was apparently removed in small quantities 
from November 1993 to June 1994. During this period 
a prototype gravity wave laser interferometer at BNO 



detected unapproved transport of materials from under- 
ground. After discovery of the theft, all of the gallium 
was cleaned, additional security controls for access to 
the gallium were instituted, and SAGE resumed oper- 
ation. As this period of time has some uncertainty with 
respect to experimental control, it is singled out for sep- 
arate treatment, and is not included in our best estimate 
for the neutrino capture rate. 

An experiment using a 517 kGi ^^Cr neutrino source 
pO| , pl| began in late December of 1994 and continued 
until May 1995. We refer to all data after January 1995 
as SAGE III, with a special designation of SAGE III.l 
for solar neutrino extractions during the Gr experiment. 

Table g summarizes the data period designations. The 
exposure times and other data for all runs of SAGE that 
are potentially useful for solar neutrino capture rate de- 



termination are given in Table III 



III. EXTRACTION OF Ge FROM Ga 

The extraction and concentration of germanium in the 
SAGE experiment consists of the following steps. 

1 . Ge is extracted from the Ga metal into an aqueous 
solution by an oxidation reaction. 

2. The aqueous solution is concentrated. 

(a) Vacuum evaporation reduces the volume of 
aqueous solution by a factor of 8. 

(b) Ge is swept from this solution as volatile 
GeCU by a gas flow and trapped in 1 1 of de- 
ionized water. 

(c) A solvent extraction is made from the water 
which concentrates the Ge into a volume of 
100 ml. 

3. The gas GeH4 is synthesized, purified, and put into 
a proportional counter. 

The average extraction efficiency from the Ga metal to 
GeH4 was 77% before 1997 and 87% thereafter. Each of 
these steps will now be briefly described and this section 
concludes with a description of the evidence that the ex- 
traction procedure does indeed remove germanium with 
high efficiency. 



A. Chemical extraction procedure 

1. Extraction of Ge from metal Ga 

A procedure for the extraction of Ge from metallic Ga 
was first investigated at Brookhaven National Labora- 
tory |2^. It is based on the selective oxidation of Ge in 
liquid Ga metal by a weakly acidic H2O2 solution. This 
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TABLE III. Parameters for all 67 runs that arc potentially useful for solar neutrino measurement. The efficiency values 
include the reduction for the energy cut (usually 0.9815) and for the rise time cut (usually 0.95). L-peak efficiencies are only 
given for those extractions for which wave form data were available. Peak ratio values, whose uncertainty is approximately 
±0.02, are only listed for those extractions for which '^'^^Cd calibrations were made. This ratio is 1.0 for unpolymerized counters. 



Mean Exposure Ga 



Exposure 


exposure 






Extraction 


Counter 


Pressure 


Percent 


O p erating 


Counting 


- p c ak 




Peak 






( H avs 1 


(tons) 


efficiency 




(mm fTcfl 


GeH4 


voltage 


system 


efficiency 


efficiency 




Tan Qf) 


1 990 040 


42.0 


28.67 


0.78 


Ni 1 


604 


28.0 


1230 


2 


0.333 






Feb. 90 


1 990 1 "^9 


30.0 


28.59 


0.79 


LA12 


635 


53.0 


1450 


5 


0.249 










26.0 


28.51 


0.81 


Ni 1 


640 


25.0 


1238 


2 


0.343 






Apr. 90 




19.0 


28.40 


0.76 


LA24 


850 


30.0 


1430 


5 


0.335 






Tn]v QO 


1 999 '^4^ 


21 


21 01 


78 


Ni 1 


524 


19 3 


1130 


2 


327 








1991 4^^'^ 


53 


27 43 


82 


LA74 


715 


28 


1300 


2 


334 






July 91 


-Lyy ± . ooy 


23 


27 37 


66 


LA77 


710 


24 


1300 




320 








1991.622 


26.3 


49.33 


0.78 


RD2 


570 


34.0 


1700 


5 


0.250 






Sep 91 


1991 707 


27.0 


56.55 


0.78 


LA40 


935 


40.0 


1630 


2 


0.338 






Nov. 91 


1QQ1 #*79 


26.0 


56.32 


0.81 


LA46 


108 


30.0 


1746 


3 


0.339 






Dec. 91 


1QQ1 Q4J^ 


26.8 


56.24 


0.79 


LA51 


870 


27.0 


1394 


2 


0.336 






Feb. 92-1 


1 QQ2 1 3S 


24.5 


43.03 


0.80 


LA71 


666 


12.0 


1110 


2 


0.322 






Feb. 92-2 


1 9Q9 1 

XJ7C7^ . XOO 


24.5 


13.04 


0.80 


LA50 


640 


30.0 


1165 


2 


0.305 






iVl ObL , f ^ 


1 999 914 
X y y^ . ^ xt: 


20 9 


55 96 


78 


LA46 


810 


20.5 


1292 


2 


316 






Apr 92 


1 999 9S4 


23 5 


55 85 


83 


LA51 


815 


23.0 


1386 


2 


333 






iVl dV ^ 


1 999 

X f C7 . O O 


27 5 


55 72 


67 


LA95 


675 


69 


1620 


2 


282 






Qpri 99 


1 999 700 


116 8 


55 60 


53 


LAllO 


720 


21 


1311 


3 


338 


0.322 




Oct 92 


1 999 790 


27 2 


55 48 


0.83 


LAlll 


725 


25.0 


1391 


3 


0.341 


0.327 




Nov. 92 


1QQ9 R71 


26.7 


55.38 


0.81 


LA105 


730 


23.0 


1351 


3 


0.315 


0.297 




Dec. 92 




24.3 


55.26 


0.85 


LA116 


740 


26.0 


1406 


3 


0.325 


0.315 


1.04 


Tan Q*^ 

J dill, yo 


1 QQ*^ O^Q 


32 3 


55 14 


76 


LAllO 


770 


25 


1412 


3 


342 


314 




Feb. 93 


Xi7^i_J. X X^J 


23.0 


55.03 


0.79 


LA107 


730 


24.0 


1336 


6 


0.315 






Apr. 93 


1 993 2S1 

X . if O X 


26.6 


48.22 


0.83 


LAlll* 


710 


23.0 


1352 


3 


0.322 






iV± Oby C/ O 


1 99*^ '^fi4 


30 9 


48 17 


82 


LA116 


705 


16 


1210 


3 


327 




1.04 


June 93 


1993.454 


30.4 


54.66 


0.80 


LAllO 


740 


24.0 


1352 


3 


0.338 


0.313 




Tnlv 9S 


1 99S 5S7 


27.9 


40.44 


0.80 


LAlll 


675 


22.0 


1266 


3 


0.353 






Alio- 9"^ 1 


1 99'^ fi'il 
xy c7o . uox 


34 


40 36 


79 


LA107 


680 


12.0 


1210 


Q 


317 




1.00 




1 99*^ fi9^{ 

X^yO.VjirO 


63.8 


14.09 


0.51 


A9 


765 


12.0 


1130 


Q 


0.322 




1.20 




1QQ3 74Q 


13.0 


14.06 


0.79 


A12 


750 


14.0 


1224 


6 


0.333 




1.00 


Oct. 93-2 


1 993 ROO 


34.7 


14.10 


0.80 


LAlll* 


710 


15.0 


1162 


3 


0.328 


0.309 


1.03 


Ort 9'^-'^ 

WLL. yo o 


1 99*^ SI 9 

XC7C70. O X^ 


24.6 


14.02 


0.84 


LA116 


665 


14.0 


1184 


3 


323 


0.299 


1.04 


Nov. 93-1 


1 99"^ 


14.0 


14.07 


0.87 


LA119 


665 


13.0 


1113 


3 


0.321 


0.316 


1.08 


Nov 9*^ 9 


1 99*^ S44 


53 4 


26 16 


52 


LAllO 


675 


9 


1094 


3 


340 


0.326 


1 00 


Dec. 93-1 


1 99S 9S(S 


30.5 


26.13 


0.78 


A19 


760 


12.0 


1287 


3 


0.336 




1.08 


Dec. 93-2 


1 99S 9S9 


39.9 


28.05 


0.80 


LAlll 


690 


12.0 


1230 


3 


0.345 


0.331 


1.02 


Jan. 94-1 


1 994 048 


42.2 


26.67 


0.82 


LA107 


760 


12.0 


1196 




0.328 




1.00 


Jan. 94-2 


1994 0.^1 


41.1 


27.44 


0.80 


LAlll* 


750 


12.5 


1065 


3 


0.308 




1.04 


Feb. 94 


1 994 1 ^7 

X3C7'±. XO 1 


28.0 


54.01 


0.64 


LA116 


600 


15.0 


1090 


3 


0.312 


0.326 


1.04 


Mar. 94 


1994.218 


31.0 


53.94 


0.78 


LA105 


625 


10.0 


1190 


3 


0.309 


0.311 


1.00 


Apr. 94 


1994 2S3 


22.5 


53.88 


0.73 


LAllO 


685 


27.0 


1331 


3 


0.328 


0.335 


1.00 


1 VI d y Ty ■4:— O 


X f CJ'X. O 1 'X 


32 9 


26 99 


85 


LAlll 


610 


17.0 


1215 


3 


329 


343 


1 00 


Tnlv 94 


1 qQ4 f^f^l 

X t7 . U X 


31 3 


50 60 


80 


LA107 


620 


22 


1236 


3 


301 


269 


1 00 


Alio- Qzl 


1 994 fi'i4 


31 


50 55 


80 


LA105 


655 


13 


1196 


3 


312 


0.307 


1 05 


Ron 94-1 


1994.722 


33.2 


37.21 


0.76 


A13 


695 


18.0 


1270 


3 


0.334 


0.319 


1.07 


Oct 94 


1 QQ4 7QQ 


28 8 


50 45 


76 


Al9 


695 


25 


1375 


3 


334 




1 06 




1 994 


31.0 


50.40 


0.79 


LA113 


685 


28.5 


1383 


3 


0.306 


0.314 


1.05 


Dec. 94 


1 994 9^1 

xyyi. fox 


21.0 


13.14 


0.80 


A12* 


610 


16.5 


1184 


g 


0.310 




1.02 


Mar. 95 


1995.209 


42.5 


24.03 


0.92 


A28 


690 


18.5 


1222 


g 


0.321 




1.00 


July 95 


1995.538 


19.9 


50.06 


0.86 


LA107 


635 


30.0 


1333 


3 


0.298 


0.317 


1.01 


Aug 95 


1 99 


46 7 


50 00 


70 


Al2 


710 


17 


1260 


3 


325 


312 


1 01 


Sep. 95 


1 995 749 


28.8 


49.95 


0.67 


LA46 


645 


37.0 


1382 


3 


0.283 


0.294 


1.02 


Oct. 95 


1995.807 


18.7 


49.83 


0.49 


A19 


680 


18.5 


1248 


3 


0.319 


0.294 


1.08 


Nov. 95 


1995.875 


25.8 


49.76 


0.89 


A9 


685 


33.0 


1429 


3 


0.310 


0.294 


1.21 


Dec. 95-2 


1995.962 


32.7 


41.47 


0.73 


LA113 


725 


18.5 


1271 


3 


0.319 


0.278 


1.00 


Jan. 96 


1996.045 


29.7 


49.64 


0.77 


A12 


715 


24.0 


1340 


3 


0.321 


0.310 


1.00 


May 96 


1996.347 


49.9 


49.47 


0.75 


LA116 


685 


21.5 


1295 


3 


0.320 


0.319 


1.00 


Aug. 96 


1996.615 


45.0 


49.26 


0.77 


A13 


675 


23.0 


1332 


3 


0.327 


0.330 


1.09 


Oct. 96 


1996.749 


45.8 


49.15 


0.83 


LA116 


635 


15.0 


1185 


3 


0.318 


0.319 


1.02 


Nov. 96 


1996.882 


48.7 


49.09 


0.78 


A12 


720 


21.5 


1308 


3 


0.323 


0.306 


1.00 


Jan. 97 


1997.019 


49.8 


49.04 


0.85 


LA113 


700 


29.0 


1372 


3 


0.308 


0.295 


1.00 


Mar. 97 


1997.151 


44.9 


48.93 


0.93 


A13 


650 


23.5 


1339 


3 


0.323 


0.335 


1.08 


Apr. 97 


1997.277 


42.9 


48.83 


0.90 


LA116 


670 


29.0 


1360 


3 


0.313 


0.320 


1.02 


June 97 


1997.403 


45.6 


48.78 


0.87 


A12 


675 


24.5 


1320 


3 


0.314 


0.314 


1.00 


July 97 


1997.537 


45.9 


48.67 


0.91 


LA51 


690 


15.5 


1242 


3 


0.321 


0.312 


1.03 


Sep. 97 


1997.671 


46.4 


48.56 


0.75 


A13 


650 


25.0 


1318 


3 


0.322 


0.335 


1.04 


Oct. 97 


1997.803 


45.0 


48.45 


0.83 


LA116 


635 


23.5 


1318 


3 


0.328 


0.327 


1.03 


Dec. 97 


1997.940 


47.0 


48.34 


0.88 


A12 


710 


27.0 


1382 


3 


0.318 


0.306 


1.00 
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FIG. 1. Chemical reactor for extraction of Ge from Ga. 



method was developed and fully tested in a 7.5-ton pi- 
lot installation at the Institute for Nuclear Research . 
The final procedure extracts Ge with high efhciency and 
dissolves only a small amount of Ga p^j25[ |. 

The Ga at BNO is contained in chemical reactors, each 
of which is able to extract from as much as 8 tons of Ga. 
The reactor (Fig. ^ is a 2-m'^ Teflon tank with ~ 40- 
mm-thick walls to which band heaters are attached. The 
Teflon tank is placed inside a secondary stainless steel 
tank. The Ga can be stirred with a motor that can turn 
an internal mixer at up to 80 rpm. A specially designed 
set of vanes are attached to the inside cover of the reactor 
that serve to completely disperse the extraction reagents 
(density 1.0 kg/1) throughout the liquid Ga (density 6.1 
kg/1). The vanes are made from Teflon and the stir- 
rer and cover are Teflon lined. A glass viewport in the 
reactor cover enables one to see the extremely vigorous 
mixing action. Ten such reactors are installed at BNO 
which are connected with a system of heated Teflon tub- 
ing and a Teflon pump that can transfer Ga between 
reactors. A system of glass- Teflon dosing pumps can put 
a measured volume of reagents into any reactor, and a 
vacuum suction device made from Teflon, glass, and zir- 
conium extends to the Ga surface to remove the reagents. 
The filling of a reactor with reagents and the stirring are 
controlled by an automated system. 

Each measurement of the solar neutrino flux begins 
by adding to the Ga approximately 700 ^g of stable Ge 
carrier (distributed equally among all of the reactors) in 
the form of a solid Ga-Ge alloy with known Ge content 
2 X lO^"' mass %). The reactor contents are then 
stirred so as to thoroughly disperse the carrier through- 



out the Ga metal. After a typical exposure interval of 
4-6 weeks, the Ge carrier and any additional Ge atoms 
produced by solar neutrinos or other processes are chem- 
ically extracted from the Ga. 

The efficiency of Ge extraction depends on a number 
of parameters. Since the efficiency falls rapidly as the Ga 
temperature increases, we begin to extract with the Ga 
at 30.0-30.5 °C, just slightly above its freezing tempera- 
ture (29.8 °C). The efficiency increases with an increase 
in the amount of oxidizing agent (H2O2), but this has the 
detrimental effect of dissolving more Ga. The efficiency 
also depends on the volume of aqueous phase which de- 
fines the time of later concentration of Ge, the most time 
consuming part of the entire extraction process. Taking 
into account all of these factors, a procedure was devel- 
oped which extracts about 85% of the Ge and dissolves 
only 0.1% of the Ga. 

The extraction solution for a reactor containing 7.5 
tons of Ga consists of 200 1 of de-ionized water, 5 1 of 
7 M HC1,| and 16 1 of a 30% solution of H2O2. All com- 
ponents of this solution are purified so their Ge content is 
negligible. Immediately after the reagents are added, re- 
actor stirring starts at a speed of 70 rpm. As the mixture 
is intensively stirred, the gallium turns into fine droplets 
which are covered with a Ga oxide film. This film pre- 
vents fusion of the droplets and holds the Ga as an emul- 
sion ||2^,^. The dissolved Ge in the Ga migrates to the 
surface of the droplets where it is oxidized and incorpo- 
rated into the oxide film. Because of the highly exother- 
mic oxidation reaction, the Ga temperature rapidly rises. 
After approximately 25 min, the II2O2 has been con- 
sumed; the Ga temperature plateaus and the emulsion 
spontaneously breaks down. To dissolve the oxide con- 
taining Ge, the extraction procedure is finished by adding 
45 1 of 7 M HCl (cooled to -15 °C) and stirring for 1-2 
min. The Ga temperature at the end of this extraction 
process is ~ 50 °C. 

The extraction solution is immediately decanted and 
sent to the first step of concentration, which is evapora- 
tion. The Ga in each reactor is then washed by adding 
20 1 of 0.5 M HCl. This solution is stirred with the hquid 
Ga for about 1 min, is decanted out, and is added to the 
previous extraction solution. Finally, to prevent oxida- 
tion of the Ga during the interval between extractions, 
a solution of 0.5 M HCl is added to the reactor and left 
there until the next extraction. 



2. Vacuum evaporation of extraction solutions 

Extraction is made sequentially from one reactor to the 
next. All the extraction solutions, whose total volume is 



*The symbol M stands for the amount of substance concen- 
tration in moles per liter. 
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2200 1 for 60 tons of gallium, arc combined at the evap- 
oration step, which is carried out in a glass recirculation 
apparatus with a steam- heated active volume of 70 1. As 
the evaporation proceeds, the acidity of the evaporated 
solution increases. Ge is volatile from concentrated chlo- 
ride solutions, so the evaporation is stopped when the 
volume of solution reaches 250-270 1, before loss of Ge 
can begin. The average time for evaporation is 15 h. 

3. Sweeping 

The next step is based on the volatility of GeCU from a 
concentrated solution of HCl. The evaporated extraction 
solution, which contains 250 g of Ga/1 in the form of 
chloride, is transferred to glass vessels with a volume of 
200 1. These vessels are part of a sealed gas flow system. 
The HCl concentration is raised to 9 M by adding purified 
12 M HCl and an air flow at 1.0 m'^/h is initiated. Ge 
is swept as GeCU from this 50 °C acid solution through 
a counter-current scrubber where the GeCl4 is absorbed 
in a 1.0 1 volume of de-ionized H2O. The amount of Ge 
remaining in the solution C{t) falls exponentially: C{t) — 
C(0) exp[— 1.84y(f)] where V is the volume of sweep gas 
in m^. The duration of sweeping is usually 2.5 h which 
gives 99% Ge extraction efficiency. At the end of the 
sweep the acidity of the absorber solution is in the range 
of 4.0 M to 4.2M, which excludes loss of Ge. 

4- Solvent extraction 

A solvent extraction is then carried out to further con- 
centrate the Ge. This procedure is based on the high 
distribution coeSicient of Ge between an acidic water so- 
lution and an organic solvent, such as CCI4. To achieve 
an optimal acidity (8.5 M), the appropriate amount of 
purified 12 M HCl is added to the solution obtained from 
sweeping. The Ge is first extracted into CCI4 and then 
is back extracted into low-tritimn H2O. This process is 
repeated 3 times. To remove the residual CCI4, a very 
small amount of hexane is added to the organic phase at 
the last step of the final back extraction. The final traces 
of hexane are removed by heating the solution at 90 °C 
for 40 min. This results in the Ge being concentrated in 
a volume of 100 ml of low- tritium H2O. 



5. Germane synthesis 

The final step of the extraction process is to synthesize 
germane (GeH4) which is used as a 20%-30% fraction of 
the counting gas in a proportional counter. NaOH is 
added to the 100-ml water solution to adjust the pH to 
the range of 8-9, and the solution is placed in a reac- 
tion flask on a high- vacuum glass apparatus. Any air 
is swept out of the solution and the connecting piping 




30 - — 
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Time after extraction start (min) 

FIG. 2. Ga temperature using two extraction procedure be- 
gun in 1997. The extraction reagents are added 1 min before 
time zero. Extraction begins when the mixer is started. The 
two vertical lines about 15 min after the start of extraction 
arc when the HCl is added to dissolve the oxide containing 
Ge. 

with a He flow and 2 g of low-tritium NaBH4 dissolved 
in 40 ml of low-tritium H2O is added. The mixture is 
then heated to 70 °C, at which temperature the Ge is 
reduced by the NaBH4 to make GCH4. The H2 gener- 
ated by the reaction and the flowing He sweep the GeH4 
onto a Chromosorb 102 gas chromatography column at 

— 196°C where it is trapped. When the reaction is fln- 
ished the column temperature is raised to — 35 °C and 
the GeH4 is eluted with He carrier gas. It is then frozen 
on another Chromosorb 102 trap at —196 °C where most 
of the He is pumped away. The GeH4 is then transferred 
with a mercury-filled Toepler pump to a glass bulb at 

— 196°C where any residual He is pumped away. The 
Toepler pump is used again to transfer the GCH4 to a 
calibrated stem, where the GeH4 volume is measured. 
During this transfer the temperature of the bulb is held 
at — 142°C so as to minimize Rn. A measured quantity 
of old low-background Xe is added and this gas mixture 
is inserted into a miniature proportional counter. The 
counter has been evacuated at 10~^ torr and baked at 
100 °C for at least 6 h. 



6. Modified procedures for SAGE III 

a. Extraction from Ga. At the beginning of 1997, the 
extraction procedure was modified to a two-step extrac- 
tion process. In the first step the volume of reagents 
added to each reactor is reduced from the values given 
previously by a factor of 2. The remaining steps in the 
removal of Ge from the Ga proceed the same as previ- 
ously described, but only require about 15 min because 
of the reduced H2O2 volume. This first step extracts 
about 75% of the Ge from the Ga, dissolves 0.05% of 
the Ga, and raises the Ga temperature to about 40 °C 
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(Fig. |[ lower curve) . After the first extraction from each 
reactor, a second extraction is carried out in the same 
order using the same volume of reagents as in the first 
extraction. By the time the second extraction begins the 
Ga has cooled to 37 °C and an additional drop of 1.5- 
2°C occurs when the new reagents are added. Since the 
initial Ga temperature is now elevated, the efficiency of 
the second extraction is less than the first, and averages 
70%. Again 0.05% of the Ga is dissolved and the final Ga 
temperature is 49 °C at the end of this second extraction 
(Fig. ^, upper curve) . This modified procedure results in 
a total efficiency of Ge removal from the Ga in excess of 
90%, but both procedures dissolve the same total amount 
of Ga (0.1%). 

b. Evaporation of extraction solutions. The vac- 
uum evaporation was modified at the beginning of 
SAGE in. Instead of stopping the distillation before the 
Ge volatilizes, the Ge is allowed to evaporate, at which 
time collection of Ge in the condensate is begun. Evapo- 
ration is continued until all the Ge has been transferred 
to the condenser. The condensate is then further evap- 
orated until its acidity is 4.5 M. This solution, whose 
volume is about 130 1, is transferred to the sweeping ap- 
paratus, 12 M hydrochloric acid is added to obtain 9 M 
acidity, and the Ge is swept out in the same manner as 
for SAGE I and II. An important advantage of this new 
method is that the solution that results from sweeping is 
pure 9 M HCl, free from Ga or Ge, so it can be used in 
later extractions. These chemical technology modifica- 
tions in SAGE III increase the efficiency of Ge extraction 
by 6%-7%, decrease the average duration of concentra- 
tion by 3-4 h, and reduce the consumption of concen- 
trated HCl by 2.5 times. 



B. Chemical extraction efficiency 



C. Tests of the extraction efficiency 

The Ga experiment relies on the ability to extract a 
few tens of atoms of ''^Ge from 5 x 10^^ atoms of Ga. 
To measure the efficiency of extraction, about 700 ^g of 
stable Ge carrier is added to the Ga at the beginning 
of each exposure, but even after this addition, the sep- 
aration factor of Ge from Ga is still 1 atom in 10"'^^. In 
such a situation one can legitimately question how well 
the extraction efficiency is known. We have performed 
auxiliary measurements to verify that this efficiency is 
well established, and briefly describe these tests in this 
section. 



1. Cr experiment 

The most direct experiment of this type involved the 
irradiation of Ga with the 747-keV neutrinos from an 
artificial source of ^^Cr [^0l|2l| ]. Eight exposures of 13 
tons of Ga were made to a 517 kCi ^^Cr source. The ^^Ge 
atoms were extracted by our usual chemical procedure 
and their number determined by counting. The ratio R 
of the measured neutrino capture cross section [P0|j21| to 



the theoretically calculated cross sections of Bahcall 
and Haxton [g8| was 



I Pll to 

5r|i]] 



R 



f^measured 



(3.1) 



^^theoretical 

0.95 ±0.12 (expt) toolv (theor) (Bahcall), 
0.87 ± 0.11 (expt) ±0.09 (theor) (Haxton). 



With either of these theoretical cross sections, R is con- 
sistent with 1.0, which implies that the extraction effi- 
ciency of ^^Ge atoms produced in Ga by the neutrinos 
from ^^Cr is the same as that of natural Ge carrier. 



The total efficiency of extraction of Ge is given by the 
ratio of the Ge content of the synthesized germane to 
the Ge present in the reactors at the beginning of the 
exposure interval. As a check, the amount of extracted 
Ge is also determined by atomic absorption analysis of a 
small fraction of the solution used in the GeH4 synthesis. 
The extraction efficiency prior to 1997 was typically 80%. 
The modified extraction procedure initiated in 1997 gives 
about a 10% higher overall efficiency. The extraction 



efficiency for each run is given in Table HI. 

Since each extraction leaves 10%-20% of the carrier Ge 
still present in the Ga, it is customary to make a second 
extraction within a few days after the first. This second 
extraction removes most of the residual Ge so that the 
Ge content of each reactor is well known after the carrier 
Ge is added. Occasionally a third extraction is made to 
totally deplete the Ge content. The extracts from these 
additional extractions are usually processed in the same 
manner as for the solar neutrino extraction, including 
counting of the synthesized GeH4. 



2. Ga(n,'y) experiment 

To test the possibility that atomic excitations might tie 
up ^"'^Ge in a chemical form from which it would not be 
efficiently extracted, the radioactive isotopes ™Ga and 
''■^Ga, which beta decay to ^"Ge and ""^Ge, were pro- 
duced in liquid gallium by neutron irradiation. The Ge 
isotopes were extracted from the Ga using our standard 
procedure. The number of Ge atoms produced was deter- 
mined by mass spectroscopic measurements [ p^ and was 
found to be consistent with the number expected based 
on the known neutron flux and capture cross section, thus 
suggesting that chemical traps are not present. 



3. Removal of Ge 

Further evidence that the extraction efficiency is well 
understood came from monitoring the initial removal 
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from the Ga of cosmogenically produced ^^Ge. This nu- 
chde was generated in the Ga as it resided outside the 
laboratory exposed to cosmic rays. When the Ga was 
brought underground, the reduction in the ^^Ge content 
in the initial extractions was the same as for the Ge car- 
rier. 



Gas fill tube 




FIG. 3. Schematic view of a proportional counter. 



4- Ge carrier 

A special Ge carrier was produced which contained a 
known number of ^^Ge atoms. This carrier was added to 
a reactor holding 7 tons of Ga, three successive extrac- 
tions were carried out, and the number of ''^Ge atoms in 
each extraction was determined by counting. The results 
verified that the extraction efficiencies of the natural 
Ge carrier and ^^Ge track each other very closely. 

IV. COUNTING OF ''^Ge 
A. General overview 

Once the ^^Ge is isolated internally in the proportional 
counter, its decay must be identified. ''-^Ge decays solely 
by electron capture to the ground state of ''^Ga with a 
half-life of 11.43 days @. The probabilities of K, L, 
and M capture are 88%, 10.3%, and 1.7%, respectively 
K capture gives Auger electrons with an energy of 
10.367 keV (41.5% of all decays), 9.2-keV x rays accom- 
panied by 1.2-keV Auger electrons from the subsequent 
M-L transition (41.2% of all decays), and 10.26-keV x 
rays accompanied by 0.12-keV Auger electrons (5.3% of 
all decays). L and M capture give essentially only Auger 
electrons with energies of 1.2 keV and 0.12 keV, respec- 
tively 1^ . The proportional counter observes the Auger 
electrons and, with considerably less efficiency, the x rays 
emitted during the relaxation of the atomic electron shell. 
As a result, about the same fraction of events occur in 
the L and K peaks. 

These low-energy Auger electrons and x rays produce 
a nearly pointlike ionization in the counter gas. This 
ionization will arrive at the anode wire of the propor- 
tional counter as a unit, resulting in a fast rise time for 
the pulse. In contrast, although a typical /3 particle pro- 
duced by a background process may also lose 1-15 keV 
in the counter gas, it will leave an extended trail of ion- 
ization. This ionization will arrive at the anode wire 
distributed in time according to its radial extent in the 
counter, which usually gives a pulse with a slower rise 
time than for a ^^Ge event. The identification of true 
^"'^Ge events and the rejection of background events are 
thus greatly facilitated by using a two parameter analy- 
sis: a candidate ^-"^Ge event must not only fall within the 
appropriate energy region, but must also have a rise time 
consistent with pointlike ionization. 



To properly determine the background rate it is neces- 
sary to count each sample for a long time after any ^^Ge 
has decayed. We endeavor to begin to count as soon as 
possible after extraction and to continue counting for at 
least 160 days. Since the number of high-quality low- 
background counters and of available counting channels 
is limited, runs are occasionally ended before the desired 
counting duration is met to permit another run to begin. 
Further, since many counters are measured in a common 
system, counting time is frequently lost for calibration or 
for counter installation or removal. 

This section continues with a discussion of how the 
proportional counters are made, how their counting effi- 
ciency is determined, and how they are calibrated, and 
concludes with a description of the counting electronics. 

B. Proportional counters 

The design and construction of the proportional coun- 
ters are based on the experience gained in the CI exper- 
iment. They are made only from materials that are ra- 
dioactively clean, are assembled in a clean environment, 
are only exposed to high levels of radioactivity during ef- 
ficiency measurement, and are always counted for back- 
ground before use in a solar neutrino extraction. 

Although several different types of counters were used 
at the beginning of the experiment, all counters used 
since the extraction of September 1991 are of a common 
type, shown in Fig. |[ The counter bodies are fabricated 
by a glassblower from Heraeus Amersil transparent syn- 
thetic fused silica (Suprasil). The main body is 10 cm 
long with an 8-mm outer diameter. One end is open for 
insertion of the cathode but can be sealed with a flared 
plug. Three tubes are attached to the other end — one 
tube with 2-mm inner diameter is used for insertion of 
the filling gas; the other two tubes are capillaries for the 
cathode and anode electrical feedthroughs. A 2-mm hole 
is made in the counter body near its center over which a 
very thin piece of blown silica is sealed. There is a cor- 
responding hole in the cathode at this position so that x 
rays from external sources can pass through this window 
into the counter gas for calibration. 

The main body of the counter is large enough to hold 
a snug-fitting zone-refined iron cathode sleeve, whose di- 
mensions are approximately 5 mm diameter, 5 cm length, 
and 1/3 mm thickness. The cathodes are individually 
machined to fit each counter body, making sure that there 
is sufficient space between the cathode and the body to 
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permit the counter to be heated to at least 100 °C for 
bakeout of impurities. The iron is drilled and cut to 
length using only new tools and ultrapure hexane as lu- 
bricant. 

The major component of the counters, Suprasil, has 
a total metallic impurity content of < 1 ppm by weight 
and OH and equivalent H2O contents of ^ 10"^ ppm. The 
cathode material typically has less than 1 ppm metallic 
impurities, except for copper which is present at ~ 7 
ppm. 

The first step in counter fabrication is a helium leak 
test of the seals and thin calibration window of the 
counter body. All parts of the counter are then thor- 
oughly cleaned: the silica parts are soaked overnight in 
aqua regia, etched briefly in hydrofluoric acid, thoroughly 
washed in high-purity water, and dried in an oven at 
slightly above 100 °C. The cathodes are washed in hexane 
in an ultrasonic bath, baked, and dried under vacuum for 
approximately 24 h at 500 °C. After cleaning, all counter 
parts are handled only with gloves and clean tools. 

The final steps of counter fabrication take place inside 
a laminar flow clean bench. Under a microscope, a 25- 
/im wire of high-purity tungsten is spot welded to the 
cathode and then threaded through a thin capillary to 
the outside of the counter, where an external lead pin is 
connected. Again under a microscope, a 12.5-/im tung- 
sten anode wire is threaded through the second capillary, 
through the center of the cathode sleeve, and welded to a 
50-/im tungsten spring wire held in place at the end of the 
counter by the Suprasil end plug. With the anode and 
cathode wires held taut in the capillaries, the electrical 
connections to the external leads are made with a small 
dab of conducting epoxy injected into the end of the cap- 
illary with a hypodermic needle. With the wires still held 
taut, the quartz end plug is gently welded in place by a 
glassblower, and then (with the counter filled with '^0.1 
atm hydrogen to prevent oxidation of the thin wires), 
the capillaries are heated and sealed around the cathode 
and anode wires. The counters are then tested for gas 
tightness, evacuated and baked for > 72 h, purged, and 
filled for testing with P-10 counter gas (90% argon, 10% 
methane) . 

Counters are tested at the time of fabrication for sta- 
bility, gain, and resolution. Counter background rates 
are measured at Baksan and are in the range of 0.1/day 
(0.07/day) in the ^^Ge L-peak (K-peak) candidate re- 
gions. 

C. Measurement of proportional counter efficiency 

This section gives a general description of the meth- 
ods of counter efficiency measurement, shows how these 
methods are applied to determine the L- and K-peak ef- 
ficiencies of several typical counters, and presents how 
the counting efficiency of the solar neutrino extractions 
is determined. 



1. Measurement methods 

Two different techniques and three different isotopes 
are employed: ^^Ar to measure volume efficiency, and 
^^Ge and ^^Ge to measure the L- and if-peak efficiencies. 

The first method uses ^^Ar to measure the volume effi- 
ciency, which we define as the probability that the decay 
of a radioactive atom in the gas phase in the volume of 
a proportional counter will produce a detectable pulse. 
The ■^''Ar source is produced by the (n, a) reaction on 
^"Ca using fast neutrons from the research breeder reac- 
tor of the Institute of Physics and Power Engineering in 
Obninsk. The extracted '^^Ar is purified on a Ti getter 
and then mixed with 90% Ar plus 10% CH4. A small 
sample of this mixture is placed into the counter under 
test, the counter high voltage is set so that the '^^Ar L 
peak is at least one-quarter scale on the energy analog- 
to-digital converter (ADC), and an energy spectrum is 
measured. The gas sample is then transferred with very 
high efficiency ^'transfer (> 99.5%) to a counter that was 
specially constructed for these measurements. It is 20 
cm in length with an internal diameter of 4 mm. It has 
a deposited carbon film cathode, shaped ends to mini- 
mize end effects, a volume of 2.5 cm'^, and a volume effi- 
ciency of (99.5±0.2)%. Additional Ar-10% CH4 is added 
to bring the pressure in this standardization counter to 
about the same value as in the test counter, and another 
energy spectrum is measured under similar conditions to 
that of the test counter. To find the position of the L 
peak, these two spectra are fit to a Gaussian plus a con- 
stant background. An energy threshold is then set at 
one-third of the peak value, equivalent to about 80 eV, 
and the total number of counts above this threshold de- 
termined by summation. This gives the count rates in 
the counter under test, i?test, and in the standardiza- 
tion counter, -Rstandard • After making minor corrections 
for background rates, the volume efficiency of the test 

counter is given by = 0.995i?tcst£'transfoi-D/^standard 

where D is the decay factor of the ^^Ar between the times 
of measurement of the two spectra. Because of the high 
and well-known transfer efficiency and standardization 
counter efficiency, the total estimated uncertainty in the 
volume efficiency of the test counter using this method is 
only 0.005, or approximately 0.6%. 

The second counter efficiency measurement method 
uses ^^Ge. A brief description is given here; for more 
details see Ref. |Q. The ^^Ge source is made by the 
{p, n) reaction on 99%-enriched ^^Ga with 7-MeV pro- 
tons from the cyclotron of the Nuclear Physics Institute 
of Moscow State University. ^^Ge is extracted from the 
gallium target, synthesized into ^^Gell4, and added to 
a normal GeH4-Xe counter filling. ^^Ge decays both by 
electron capture (64%) and by positron emission (36%). 
About 40% of the electron capture decays go to an ex- 
cited state of the daughter ^^Ga which emits a coinci- 
dent 1106-keV gamma ray. The measurements are made 
by placing a proportional counter with a ^^Ge filling on 
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TABLE IV. Volume efRciency measurements using the 
three techniques described in the text. 



2. Application and test 



Counter 



Volume efficiency measured with 



'Ar 



'Ge 



'Ge 



LA51 


0.8S7 


-|- 


0.005 






LASS 


0.876 


-|- 


0.005 


n 854- + 0015 


87Q + 01 5 


LA105 




-j- 


n OfTi 

u.wo 






LA107 


0.S74 


-|- 


0.005 






LAllO 


0.933 


± 


0.005 






LAIU 


0.94S 


± 


0.005 






LAlll* 


0.S97 


± 


0.005 


0.895 ± 0.015 


0.908 ± 0.015 


LA113 


0.875 


± 


0.005 






LA114 


0.892 


± 


0.005 


0.918 ± 0.015 


0.913 ± 0.015 


LA116 


0.901 


± 


0.005 






AS 


0.868 


± 


0.005 


0.867 ± 0.015 




A13 


0.928 


± 


0.005 






A2S 


0.893 


± 


0.005 






A31 


0.872 


± 


0.005 







Average 0.894 ± 0.025 



the axis of and 10 cm to 12 cm distant from a large Ge 
semiconductor detector which observes the gamma rays. 
Energy spectra are taken of the events produced by elec- 
tron capture decays of ^^Ge by gating the signal from the 
proportional counter with the output of a single channel 
analyzer set on the 1106-keV gamma ray. The K, L, and 
volume efficiencies are defined as the ratio of the number 
of counts in the K peak, L peak, and total spectrum, re- 
spectively, to the number of 1106-keV gammas detected 
by the large germanium detector. In these calculations, 
small corrections are made to the raw number of observed 
events because of random coincidences and background 
in the Ge detector. The uncertainty in this measurement 
method is mainly from the partial detection of Af-peak 
events. The M peak in Ge is at ~ 120 eV, a higher en- 
ergy than in Ar, but it contains a much larger fraction of 
the total number of decays (7% compared to 1.4% in Ar). 
Even though part of the M peak is detected in these Ge 
spectra, a substantial correction for the missing fraction 
of events below threshold energy is still required. The es- 
timated uncertainty in the peak efficiency is thus slightly 
less than 2.5% (or 0.008 in absolute efficiency) and 1.7% 
(0.015 in absolute efficiency) in the volume efficiency. 

The final measurement method uses ''^Ge produced by 
neutron irradiation of ™Ge. After extraction and pu- 
rification of the Ge, ^^GeH4 is synthesized, and mixed 
with Xe-GeH4. Measurements of the volume efficiency 
are then made using a similar technique to that described 
for ^^Ar. In addition, the L- and K-peak efficiencies are 
determined by integration over the peaks. The uncer- 
tainty in this measurement method is about the same as 
for the ^^Ge method. 



Table IV gives the measured volume efficiencies for 14 
counters using the measurement methods based on "^^Ar, 
^^Ge, and ^^Ge. For those counters that were measured 
with more than one isotope, the agreement is very good 
and distributed in the expected statistical manner. 

The efficiencies in the L and K peaks for four coun- 
ters measured with the ^^Ge coincidence method and for 
three counters measured with Ge are given in Table 0. 
Because different gas compositions and pressures were 
used in these counter fillings, these measurements can 
only be compared if one has a procedure for correcting 
the efficiency for the gas filling. Since our solar neutrino 
runs also have different counter fillings, such a correction 
procedure is also essential for determining the counting 
efficiency for normal extractions. 

The ""^Ge counting efficiency e{P, G), before the appli- 
cation of energy or rise time cuts, can be written in the 
general form 



e(P,G)-e,(l-/z,)i?(P,G), 



(4.1) 



where is the volume efficiency, G is the fraction of the 
counting gas that is GeH4, P is the total counter pressure 
in standard atmospheres, and fo is the fraction of peak 
events that lie outside the ±1 full width at half maximum 
(FWHM) energy window, determined empirically for our 
counters from ^^Ge and ^^Ge spectra to be 0.063 for the 
L peak and 0.202 for the K peak. Monte Carlo simu- 
lations, based on our standard counter geometry, were 
made to determine the dependence of the efficiency on 
P and G [Q. Fits to these calculations with a polyno- 
mial function give E{P, G) = A{G) + B{G)P + C{G)P'^, 
where A{G) = Aq + AiG, B{G) = Bo + ^iG, and 
C(G) = Go + GiG. This equation apphes to both the 
L and K peaks with different constants in the expres- 
sions for A, B, and G. For the L peak the constants are 
Ao = 51.0, Ai = 5.51, Bo = -15.7, Bi = 1.58, Go = 3.0, 
and Gi = 0.000113, and for the K peak the constants 
are Aq = 29.7, Ai = -8.27, Bq = 28.4, Bi = -5.02, 
Go = —6.22, and Gi = 2.27. Over the range of counter 
fillings for usual extractions, the estimated uncertainty 
in E from the Monte Carlo calculations is ±1%. 

With the aid of this efficiency formula it is now possi- 
ble to compare the measurements in Table ^ Calculated 
efficiencies for these counters in the L and K peaks, us- 
ing the volume efficiency measured with '^''Ar and Eq. 
( [4.lD , are given in columns 6 and 8 of Table The total 
uncertainty in the calculated efficiencies is estimated to 
be 1.5%, consisting of 0.6% from uncertainty in ey, 1.0% 
from uncertainty in fo, and 1.0% from the uncertainty 
in the Monte Carlo simulations. The calculated efficien- 
cies agree with the values measured with ^^Ge and ^^Ge 
within the errors of calculation and measurement. 
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TABLE V. Comparison of measured counting efficiency in a 2 FWHM wide energy window centered on the L and K 
peaks with the efficiency calculated from the efficiency formula, Eq. (|l.l|), using volume efficiency measured with ^'^Ar. The 
uncertainty in the measured efficiency is ± 0.008 and the uncertainty in the calculated efficiency is estimated to be ± 1.5%. 









GeH4 




Counting 


efficiency 




Counter 


Isotope 


Pressure 


fraction 


in 


K peak 


in 


L peak 


name 


used 


(mm Hg) 


(volume %) 


Measured 


Calculated 


Measured 


Calculated 


LA88 


^«Ge 


640 


10.6 


0.326 


0.329 


0.313 


0.327 


LA88 


"Ge 


735 


9.5 


0.332 


0.345 


0.322 


0.316 


LAIU* 


««Ge 


710 


15 


0.347 


0.345 


0.334 


0.330 


LAIU* 


"Ge 


735 


9.5 


0.358 


0.353 


0.321 


0.324 


LA114 


s^Ge 


745 


8 


0.355 


0.354 


0.316 


0.320 


LA114 


^^Ge 


909 


17.4 


0.379 


0.369 


0.320 


0.310 


A8 


««Ge 


800 


12 


0.348 


0.349 


0.310 


0.308 



3. Counting efficiency for solar neutrino extractions 

The counters used during the course of the experiment 
are listed in Table III . The counter type used for the ma- 
jority of extractions is indicated by the designation "LA" 
or "A". The second type was used for three extractions 
during 1990 and is indicated by "Ni"; the final type was 
used only for the August 1991 extraction and is indicated 
by the designation "RD" . 

The counting efficiency used for each extraction is cal- 
culated by Eq. (4.1) and is given in Table [IL The vol- 
ume efficiency of most counters has been directly mea- 
sured with '^^Ar; if a counter's volume efficiency has not 
been measured, it is assumed to equal the average of all 
measured counters. Because the analysis reported in this 
section resulted in new counter efficiencies for SAGE I, 
these revised efficiencies are given in this table and are 
used in any combined fits which include SAGE 1 data. 



D. Counter calibration 

Immediately after filling counters are calibrated 
through their side window with the 5.9-keV x rays from 
an ^^Fe source. They are recalibrated with ^^Fe after 
about 3 days of operation, and then again approximately 
every 2 weeks until counting ends. This usually gives 
more than ten ^^Fe calibrations, with at least four during 
the first month of counting while the ^^Ge is decaying. 
In addition, beginning with SAGE II, calibrations are 
usually made with a ^°^Cd source whenever an ^^Fe cali- 
bration is done. The 22-keV Ag x rays that follow ^°^Cd 
decay pass through the counter window and fluoresce the 
Fe cathode, giving the K x-ray peak from Fe at 6.4 keV. 
Although these x rays originate near the counter win- 
dow, they are absorbed throughout the counter volume, 
and thus give the average counter response. Beginning 
with the February 1993 extraction, a ^"^Cd-fSe source 
was periodically used. The Cd x rays fluoresce a Se tar- 
get whose L and K x rays enter the counter through its 
side window and give peaks at 1.4 keV and 11.208 keV. 



TABLE VI. Summary of the external calibration source 
x-ray energies. 



Energy (keV) Source 



Origin 



1.4 


i°«Cd+Se 


1.625 


sspe 


5.895 


55pe 


6.4 


i°9Cd 


11.208 


^o«Cd+Se 



Se L X ray through window 
Xe escape peak through window 
Mn K X ray through window 

K X ray from Fe cathode 
Se K X ray through window 



peaks from these calibration 
in Table VI, These various 



The energies of the 

sources are summarized 

calibration lines have been used to check the linearity 
of the energy and amplitude of the differentiated pulse 
(ADP) counting channels and to determine offsets. There 
are also Xe escape peaks with the ^"^Cd and ^°^Cd-l-Se 
sources, but these lines are usually weak and not useful 
for energy scale determination. 

The typical counter resolution measured with an ^^Fe 
source is in the range of 20%-23%. Scaling the resolution 
by the square root of the energy, this implies resolutions 
in the ^^Ge L and K peaks of 45%-50% and 15%-17%, 
respectively, values that are observed in ^-"^Ge- filled coun- 
ters operated at low voltage. 

Many calibrations are done on each counter. With 
each calibration a small fraction of the GeH4 molecules 
are broken into fragments which can be deposited on the 
anode wire near the counter window. This process, which 
we call "polymerization" , gradually increases the anode 
diameter, reduces the electric field, and gives a depres- 
sion of the apparent energy measured with an ^^Fe source 
or a -"^^^Cd-I-Se source. Polymerization (see, e.g., |^) oc- 
curs most readily at high count rates, so we maintain the 
rate below 10 events/s during calibration. A check for 
the presence of polymerization is made by comparing the 
peak positions of the 5.895-keV line from the ^^Fe source 
(which provides events only at the counter window) and 
the 6.4-keV line from the ^"^Cd source (which provides 
events over a much larger fraction of the counter volume). 
If the counter anode is not polymerized near the window 
and the energy channel is linear, the ratio of peak posi- 
tions will be 6.4/5.895 = 1.086. For each extraction the 
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ratio of the 6.4-keV to 5.895-keV peak pos itions averaged 
over all calibrations is given in Table [II relative to the 
unpolymerized value of 1.086. Most counters show little 
or no evidence of polymerization. For polymerized coun- 
ters the peak ratio is greater than 1.00 and is used to 
correct the energy scale derived from each ^^Fe calibra- 
tion. 



E. Linearity of counter gain 

Calibrations with the ^°^Cd-|-Se source have been used 
to check the predicted position and resolution of the '^^Ge 
K peak from an ^^Fe calibration. Since the 11.208-keV 
peak energy with the ^"^Cd-I-Se source is very close to the 
10.367-keV energy of ^^Ge iC-peak events, this method 
has the advantage that very little extrapolation of the 
peak position in energy is needed. Some departures from 
linearity are present in the region of the ^-'^Ge K peak. 

Measurements have been made as a function of GCH4 
fraction G, counter pressure P, and operating voltage V . 
The ratio of the peak positions is equal to the ratio of the 
energies (11.2/5.9) up to a critical voltage Fcrit = 10.5G+ 
0.6P -I- 588. Above this critical voltage, the location of 
the ^^Ge K peak [P/f(^^Ge)] can be inferred from the 
location of the ^^Fe peak [P(^^Fe)] using the formula 



TABLE VII. Specifications of counting systems 2, 6, and 3. 



Pk(^^Gc) _ 10.367 
P(55Fe) ~ 5.895 



[1 - (4.5G + 2.78)(1/ - Fcrit) x 10-*^] 

(4.2) 



where G is expressed in percent, P is in mm Hg, and V 
is in volts. The typical correction due to the nonlinearity 
of the gain is a reduction in the predicted ''^Ge peak 
position of 2%. 

This set of experiments also measured the resolution 
of the peaks from ^°^Cd-|-Se and from ^^Fe. Below a 
critical voltage, the ratio of the resolutions was equal to 
the expected value of -^/S. 9/11. 2, but above this voltage, 
given by V„it = 6G-hP/3-h824, the ^"^Cd+Se resolution 
was wider than predicted from the ^^Fe resolution. From 
these measurements the relationship between the ^-'^Ge 
i^T-peak resolution [RKiJ^Ge)] and the ^^Fe resolution 
[P(^^Fe)] was found to be 



Pg(^^Ge) 
P(55Fe) 



5.895 
10.367 



[1 + 1.5 X 10-3(F-l^erit)]. (4.3) 



Note that the value for T4rit for the resolution correction 
is not the same as for the gain correction. The typical 
correction for the K peak results in an increase in the 
predicted ''^Ge resolution of 15%. 

The correction to the gain and resolution predicted by 
these empirical formulas is accurate to about 30%. The 
nonlinearity in gain and resolution is only present at the 
higher energies. No corrections are required for the ^^Ge 
L peak because the critical voltages are much higher than 
for the K peak. 



Specification 


Sys. 2 


Sys. 6 


Sys. 3 


Number of channels 


7 


7 


8 


JN umber of channels with INal 


5 


6 


8 


Counter dynamic range (keV) 


0.4-13 


0.5-18 


0.3-18 


JNal dynamic range (keV) 


50-3000 


50-3000 


50-3000 


Max. counting rate (s~^) 


5 


1000 


1.5 


Nal coincidence window {^s) 


8 


4 


5.2 


Energy time constant (/xs) 


1 


NA 


NA 


ADP time constant (ns) 


10 


10-500 


10 


Nal time constant (/is) 


1 


0.5 


1 


Bandwidth, —3 dB (MHz) 


90 


45 


90 


Rise time, 10%-90% (ns) 


3.5 


8 


4 


JNoise, peafi to peafi (mV) 


< 10 


< 12 


< 10 


Dead time 








in acquisition mode (ms) 


200 


1 


600 


in calibration mode (ms) 


50 


1 


120 


ADC resolution (mV/ch) 


10 


1 


1 


Energy ofltset (ch) 











ADP offset of 4096 ch (ch) 


-45 to +25 





70-120 



F. Electronic systems 



As indicated in Table HI, SAGE has used several dif- 



ferent counting systems as the experiment progressed. 
Most runs of SAGE I were counted in what we call sys- 
tem 2. Since the fall of 1992, during SAGE II and III, 
most first extractions were counted in system 3. System 
6 measured a few first extractions, but most were from a 
low mass of Ga. The major specificatio ns of these vari- 
ous counting systems are given in Table VI] . Since this 
article focuses on SAGE II and III, we will mainly con- 
sider counting system 3 in the following. Some additional 
information concerning system 6 and previous systems is 
given in Appendix 

The counting systems reside in a specially designed air- 
conditioned room in the underground facility. To mini- 
mize the fast neutron and gamma ray flux, the walls are 
made from low-radioactivity concrete with an outer steel 
shell. The entire room is lined with sheets of 1 mm zinc- 
galvanized steel to reduce radio-frequency noise. Power 
to the counting electronics is supplied by a filtered un- 
interruptible power supply, with signal and power cables 
laid inside independent steel conduits. The data acqui- 
sition computers, which are in the counting room, are 
networked so that the systems can be monitored outside 
the underground laboratory. The counting room is kept 
locked and access is restricted to counting personnel. 

System 3 was moved to BNO and installed in the un- 
derground counting room in 1988. It can record events 
from up to eight counters which are placed inside the well 
of a Nal crystal that serves as an active shield (crystal, 
23 cm diameter by 23 cm height; well, 9 cm diameter by 
15 cm height). There are two layers of passive shielding. 
An inner layer of square tungsten rods (10 mmxlO mm) 
encloses the Nal and the photomultipliers are shielded 
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by Pb. All components are made from low-radioactivity 
materials which were assayed prior to construction by a 
low-background solid-state Ge detector. The preampli- 
fiers are mounted as close to the counters as possible, but 
are separated by a thick layer of copper. The counters 
are sealed nearly air tight inside the apparatus. Dry ni- 
trogen gas from evaporation of liquid nitrogen flows con- 
tinuously through the Nal well to remove Rn. The en- 
tire apparatus may be lowered with a hoist into an outer 
shield whose bottom and sides consist of 24-32 mm of 
copper, 210 mm of lead, and 55 mm of steel, and whose 
top has 34 mm of Cu and 250 mm of steel. The cav- 
ity between the inner and outer shields is also purged 
continuously by gas from evaporating liquid nitrogen; to 
preclude backstreaming the flow rate is kept below that 
in the inner shield. 

To minimize the length of the signal cables, the rack 
of counting electronics is immediately adjacent to the 
outer passive shield. The electronics is in a single rack 
designed to reduce rf interference. The block diagram 
of a single channel of system 3 is illustrated in Fig. ^. 
Briefly, the analog signal processing proceeds as follows: 
the proportional counter anode is directly connected to a 
charge-sensitive preamplifier. After further amplification 
the signal is split, with one channel going to the digital 
logic to determine that an event from that counter has 
occurred, and a second channel going to a 90 ns cable de- 
lay and then to a gated multiplexer. The signals from all 
8 counters are input to separate gates of this multiplexer 
and the appropriate gate is opened by the digital logic for 
whichever counter has seen an event. The multiplexed 
output is split into four channels: two go to a digital 
oscilloscope which records the counter wave form with 
8-bit resolution for 800 ns after pulse onset at two dif- 
ferent amplification ranges, one appropriate for the ^^Ge 
L peak and the other appropriate for the K peak. One 
of the other two signals goes to an integrating ADC to 
measure the total pulse energy; the second signal is dif- 
ferentiated with a time constant of 10 ns, stretched, and 
input to a peak-sensing ADC. This second ADC measures 
the amplitude of the differentiated pulse, called "ADP." 
Acquisition can be run in calibration or event acquisition 
modes. For each event in acquisition mode, the energy, 
ADP value, time of event, Nal time and energy, and the 
two digitized wave forms (high- and low-gain channels) 
are written to disk. 



V. SELECTION OF CANDIDATE ^^Ge EVENTS 

The counting data consist of a set of events for each 
of which there is a set of measured parameters, such as 
wave form, energy, Nal coincidence, etc. The first step of 
analysis is to sort through these events and apply various 
selection criteria to choose those events that may be from 
^^Ge. We will describe here the selection procedure for 
events measured in counting system 3; the procedure for 



system 6 is identical except there are no measured wave 
forms, so the energy is measured by an ADC and the 
ADP method is used for rise time determination. 



A. Standard analysis description 

The various steps to select potential ^^Ge events are 
the following: 

(1) The first step of event selection is to examine the 
event wave form and identify two specific types of events: 
those that saturate the wave form recorder and those that 
originate from high- voltage breakdown. Saturated events 
are mostly produced by alpha particles from natural ra- 
dioactivity in the counter construction materials or from 
the decay of ^^^Rn that has entered the counter during 
filling. Such events are easily identified and labeled by 
looking at the pulse amplitude at the end of the wave 
form. Saturated pulses have amplitude greater than 16 
keV and occur in an average run at a rate of approxi- 
mately 0.5/day. Since most such pulses are seen after 
any initial ^^^Rn has decayed, they are mainly from in- 
ternal counter radioactivity. Events from high-voltage 
breakdown have a characteristic wave form which rises 
very steeply and then plateaus. A true pulse from ^^Ge 
decay, in contrast, rises more slowly and after this initial 
rise, has a slow, but steady, increase in amplitude as the 
positive ions are collected. Breakdown pulses are identi- 
fied by determining the slope of the wave form between 
500 and 1000 ns after pulse digitization begins. 

(2) To minimize the concentration of Rn, the air in 
the vicinity of the counters is continuously purged with 
evaporating liquid nitrogen. Counter calibrations, how- 
ever, are done with the counter exposed to counting room 
air which contains an average of 2 pCi of Rn per liter. 
When the shield is closed and counting begins, a small 
fraction of the decays of the daughters of ^^^Rn can make 
pulses inside the counter that mimic those of ^^Ge. To re- 
move these false ''^Ge events, we delete 2.6 h of counting 
time after any opening of the passive shield, and estimate 
the background remov al efficie ncy of this time cut to be 
nearly 100%. See Sec. |VII D 2| for further details. 

(3) It is possible that the Xc-GcH4 counter filling may 
have a small admixture of ^^^Rn that enters the counter 
when it is filled. Most of the decays of Rn give slow 
pulses at an energy outside the ^^Ge peaks, but approxi- 
mately 8% of the pulses from Rn and its daughters make 
fast pulses in the K peak that are indistinguishable from 
those of ""^Ge. Since Rn has a half-life of only 3.8 days, 
these events will occur early in the counting and be falsely 
interpreted as ''^Ge events. Each ^^^Rn decay is, how- 
ever, accompanied by three a particles, which are de- 
tected with high efficiency and usually produce a satu- 
rated pulse in the counter. Since the radon decay chain 
takes on average only about 1 h from the initiating decay 
of 222Rn to reach ^lOpb with a 22-yr half-life, deleting aU 
data for a few hours around each saturated event removes 
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FIG. 4. Block diagram of one channel of the eight-channel system 3 counting electronics. Abbreviations: PC, proportional 
counter; DAQ, data acquisition. 
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TABLE VIII. Effect of cuts on the experimental live time 
and events for all runs of SAGE II and SAGE III that were 
counted in both L and K peaks (except May 1996). The 
results of the cut on each row include the effect of all cuts on 
preceding rows. Because the rise time cut varies with energy, 
no entry can be given for the "2-15 keV" column. 





Live 


Number of events 




time 


2-15 


L 


K 


Cut description 


(days) 


keV 


peak 


peak 


None 


4129 


4209 


1990 


821 


Shield open time cut 


4040 


3962 


1864 


785 


Saturated event time cut 


3862 


3641 


1733 


728 


Nal coincidence cut 


3862 


1275 


1106 


519 


Rise time cut 


3862 


NA 


408 


314 



most of these false ^^Ge events. We choose to delete from 
15 min prior to 3 h after each saturated pulse. The ef- 
ficiency of this cut in time is 95%. Further details are 



given in Sec. VII D 1 



(4) All events whose pulse is coincident with a NaT 
detector response are then eliminated. Since ''^Ge has 
no 7 rays associated with its decay, this veto reduces 
background from natural radioactivity. 

(5) The next step is to set the energy windows for the 
Ge L and K peaks. The measure of energy is the integral 
of the pulse wave form for 800 ns after pulse onset. The 
peak position for each window is based on the calibration 
with ^^Fe, with appropri ate co rrections for polymeriza- 
tion, as described in Sec. IV D| , and for nonlinearity, as 
described in Sec. IV E. If the peak position changes from 



one calibration to the next, then the energy window for 
event selection is slid linearly in time between the two cal- 
ibrations. The resolution at each peak is held constant 
and is set to be the average of the resolutions with ^^Fe 
for all counter calibrations, scale d to the L- or if-peak 
energy as described in Sec. |IVD[ (I n the rare cases that 
the resolution of the first ^^Fe calibration is larger than 
the average, the resolution of the first calibration is used 
throughout the counting.) Events are then accepted as 
candidates only if their energy is within ±1 FWHM of 
the central peak energy. 

(6) Finally, events are eliminated unless their rise time 
is in the range of what is expected for ^^Ge decays. For 
runs with wave form recording, the rise time is derived 
from a fit to the pulse shap e wit h an analytical function, 
as described below in Sec. |VB| . For those runs without 
wave form recording, the L peak is not analyzed and the 
ADP measure of rise time is used to set the acceptance 
window for i^-peak events. 

For the 30 runs of SAGE II and III that could be 
counted in both the L and K peaks, the effect on the 
live time of each successive cut and the total number of 



candidate Ge events that survive is given in Table VIII . 
(The run of May 1996 is excluded because the counter 
was slightly contaminated with residual "^^Ar which had 
been used to measure this counter's efficiency.) Figure ^ 




FIG. 5. Upper panel shows the energy rise time histogram 
of all events observed during the first 30 days after extraction 
for all runs that could be counted in both L and K peaks 
(except May 1996). The live time is 711.1 days. The expected 
location of the ^^Ge L and K peaks as predicted by the ^^Fe 
and ^'^''Cd calibrations is shown darkened. Lower panel shows 
the same histogram for all events that occurred during an 
equal live time interval at the end of counting. 



shows all events from these same runs that survive the 
first four cuts. Events that occurred early in the counting 
are shown in the upper panel and at the end of counting 
in the lower panel. 

Several runs were compromised and some were com- 
pletely lost due to operational failures. Failure of an 
electronic component made it impossible to use the L 
peak in the extractions of April 1993, May 1993, July 
1993, and October 1994. Similar problems made it im- 
possible to make a rise time cut in the K peak for the 
runs of June 1991, July 1993, October 1994, and October 
1997. These runs thus have a larger than normal number 
of events. If an electronic component fails that deterio- 
rates the rise time response and the failure occurs early 
in the counting, while the ''^Ge is decaying, our pohcy 
is to not use any rise time cut in the K peak and to re- 
ject this run in the L peak. If the failure occurs later, 
the rise time cut is retained and the interval of failure 
is removed from the data. Extractions in March 1993, 
January 1995, May 1995, and March 1996 were entirely 
lost due to counter failure. The extractions of Septem- 
ber 1993, September 1994-2, and July 1996 were lost be- 
cause either the counter stopcock failed or some other 
gas fill difficulty occurred. Electronic failures caused the 
loss of the extractions of September 1993-1, May 1994- 
2, and April 1995. Extractions in June 1995 were lost 
due to radioactive contamination of the counters with 
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isotopes that were being used at this time for counter 
efficiency measurement. Finally, we exclude several ex- 
tractions from one reactor that were systematic studies 
in preparation for the Cr source experiment. Since their 
mass was no more than 7.5 tons of gallium, less than one 
atom of '^^Ge is detected on the average in such runs in 
the combination of both the L and K peaks. Two-reactor 
extractions, however, whose mass is approximately 15 
tons, give on the average 1.5 ''^Ge events, sufficient to 
determine the solar neutrino capture rate, albeit with a 
large error [p6[. 



B. Rise time analysis techniques 



As described in Sees. II B and IV F, the data acqui- 



sition system electronics has evolved over the course of 
SAGE. The data from SAGE I relied entirely on a hard- 
ware measurement of the rise time. This ADP technique 
suffices well in studies of the if-peak counter response, 
but is not capable of adequately differentiating rare L- 
peak events from noise. 

Wherever possible for SAGE II, and throughout 
SAGE 111, we derive a parameter that characterizes the 
rise time from the wave form, and are thus able to present 
both L- and if-peak results. For those runs with only 
ADP data, the L peak cannot be analyzed and we present 
only if-peak data. All wave form data come from count- 
ing system 3. 
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FIG. 7. Background candidate event in K peak. Note the 
much slower fall when the pulse begins at ~ 200 ns than for 
the true ^'^Ge K-peak event in Fig. 



1. Waveform rise time determination - Tn 

Figure || shows typical pulses in the L and K peaks 
from a ^^Ge- filled counter as captured by the digitizing 
oscilloscope in system 3. There are 256 channels full scale 
on the y axis corresponding to 1.040 V (130 mV/div) for 
digitizer channel 1 and 0.160 V (20 mV/div) for channel 
2. The X axis has 1024 digitization points each with 1 
ns duration. The relevant features of the pulses are the 
base line from t = to roughly 120 ns, the dc offset 
that occurs when the gate opens at 120 ns, and the fast 
onset of the pulse at about 180 ns. The exact values of 
these times and offsets vary depending on the counting 
channel and the run; they even vary slightly from pulse to 
pulse within a given run. When determining the energy 
and rise time of the pulse, it is therefore necessary to 
determine accurately the onset of the pulse both in time 
and dc voltage level. 

By treating the trail of ionization in the proportional 
counter as a collection of point ionizations and integrat- 
ing over their arrival time at the anode, it can be shown 
||37| that the voltage output V of an infinite bandwidth 
preamplifier as a function of time t after pulse onset has 
the form 



V{0<t< Tn) 
Vit > T^^) 



= Vol ^^Hl + ±)-^] 

t-TN -) , 



Vol ln(l 



^01 -n^ ■ to 
-^ln(l 

J-N ^ 



1 



t+to 



(5.1) 



with V{t < 0) = 0, where Tn is the time duration 
over which the ionization arrives at the anode, to is a 
time inversely proportional to the ion mobility, and Vq 
is proportional to the total amount of ionization de- 
posited in the counter. The parameter T/v character- 
izes the rise time of the wave form. For the case of 
true point ionization, Tn should be near zero. When 
T/v is zero, the function reduces to the Wilkinson form 
v{t : Tat = 0) = Vbln(l + t/to). When Tn is large, the 
event is characteristic of extended ionization, and is most 
likely a background event from a high-energy (3 particle 
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traversing the counter. Figure ^ is an example of such a 
slow pulse in the K peak. 

Because this form for the pulse shape has a sound phys- 
ical basis and reasonable mathematical simplicity, wc fit 
every pulse t hat is not identified as saturation or break- 
down to Eq. (5^). To account for the fact that the pulse 
onset time tonset is not at time zero, we replace t by 
t — ionsct, and since the pulse begins at a finite voltage 
VoSset 1 we replace V hy V — T4ffsot ■ The fit is made from 
40 ns before the time of pulse onset to 400 ns after onset. 
Five parameters are determined by the fit: ionsetj ^Dfisetj 
Vo (a measure of the energy deposited during the event 
which is not used in analysis), (whose value of slightly 
less than 1 ns is approximately constant for all pulses), 
and Tjv (the rise time). 



2. Alternative wave form analysis methods 

Although we use fits to Tjv as our standard analysis 
technique, we also developed two alternative methods 
to discriminate pointlike ionization from extended-track 
ionization in the proportional counter pulses. These serve 
as checks on the event selection based on T/v- One tech- 
nique is based on a fast Fourier transform (FFT) of the 
digitized wave form. No specific functional form for the 
pulse is assumed and hence this method has the advan- 
tage that it is sensitive to potential alterations in the 
pulse shape. See Appendix |b| for further information con- 
cerning the FFT method. The second method of wave 
form analysis that was investigated also assumes no par- 
ticular form for the pulse. This method, called the "RST 
method," deconvolutes the observed wave form to find 
the initial ionization pattern in the counter. See Ap- 
pendix ^ for further details. 

Since these three techniques are sensitive to different 
characteristics of the wave form, their selection of events 
is, not unexpectedly, different. Nonetheless, when many 
data sets are considered in combination, their results for 
the overall production rate are in good agreement, which 
provides strong support for the validity of our wave form 
analysis procedure. 



3. Hardware rise rime measurement: ADP 

The amplitude of the differentiated pulse is propor- 
tional to the product of the original pulse amplitude and 
the inverse rise time. The quantity ADP/energy is thus 
proportional to the inverse rise time. Events due to low- 
energy Auger electrons and x rays that produce point 
ionization in the counter all have a fast rise time. Events 
with a slower rise time (small ADP) are due to back- 
ground pulses that produce extended ionization. Events 
with a very fast rise time (large ADP) are due to elec- 
tronic noise or high-voltage breakdown. 



Inherent in an ADP analysis is the uncertainty that 
arises from an imprecise knowledge of the offset for a 
given run. Nonzero offset occurs when the gate is opened 
after an event trigger. The electronic components which 
process the pulse are subject to small drifts in their off- 
sets that are functions of external parameters, such as 
temperature. These nonzero offsets contribute to the dc 
offset on which the event pulse rides. Our approach has 
been to extrapolate ADP vs energy plots from the ^^Fe 
calibrations using the 5.9-keV peak and the escape peak 
to obtain an offset for each calibration. Since the offsets 
are typically distributed in a Gaussian manner with a 
sigma of 1 or 2 channels, the average is a good approxi- 
mation when determining the i^T-peak selection window. 
For the L peak, however, uncertainties of a few channels 
lead to significant variations in event selection. Utilizing 
the digitized pulses, it is possible to eliminate this un- 
certainty by determining every offset on a pulse-by-pulse 
basis. 

A further disadvantage of the ADP method is that it is 
only responsive to the initial rise of the pulse. Occasional 
small pulses from high- voltage breakdown have rise time 
the same as for true L-peak ^^Ge pulses, but after their 
initial rise they turn flat, rather than gently rise as the 
positive ions are collected as with a real ^^Ge event. A 
breakdown event of this type is not distinguished from a 
^^Ge event by the ADP method, but is easily recognized 
by examining the recorded wave form long after pulse 
onset. 



C. Calibration of rise time response 

To determine the values oiT^ for true ''^Ge pulses, we 
have filled counters with typical gas mixtures (20% GeH4 
and 80% Xe at a pressure of 620 mm Hg) , added a trace 
of active ^^GeH4, and measured the pulses in each of the 
system 3 counting slots. All events inside 2 FWHM of 
the L and K peaks are then selected and the rise time 
Tjv of each event calculated with Eq. (5.1). The rise time 
values are arranged in ascending order and an upper rise 
time limit set such that 5% of the events are excluded. 
This leads to event selection limits on Tm of 0.0-10.0 ns in 
the L peak and 0.0-18.4 ns in the K peak. The variation 
with electronics channel and with counter filling, over 
the range of our usual gas mixtures, was measured to 
be approximately 1.2 ns. We choose to fix the event 
selection limits at the values given above, and include 
in the systematic error an uncertainty in the efhciency 
of ±1% due to channel and filling variations. A major 
advantage of using T/v is that the rise time limits are 
fixed and are the same for all extractions. The purpose 
of the calibrations with ^^Fe and other sources is solely 
to determine the energy scale. 

For those runs in which the ADP method of rise time 
discrimination is used, the limits for the ADP cut are 
determined separately for each run from the ^^Fe cali- 
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brations. Histograms of the values of ADP/energy for 
the events within 2 FWHM of the 5.9-keV energy peak 
are analyzed to determine the cut point for 1% from the 
fast region (to eliminate noise) and 4% from the slow re- 
gion (to eliminate background). All calibrations from a 
run are analyzed and the ADP window for event selection 
is slid linearly in time from one calibration to the next. 



VI. STATISTICAL ANALYSIS AND RESULTS OF 
SINGLE RUNS 

In this section we describe how the data are analyzed 
to determine the ^^Ge production rate. We then give the 
results for individual runs and for all runs in the L- and 
K-^eak regions. 



A. Single-run results 

The above selection criteria result in a group of events 
from each extraction in both the L- and if -peak regions 
which are candidate ''^Ge decays. To determine the rate 
at which ^^Ge was produced during the exposure time, 
it is assumed in each peak region that these events orig- 
inate from two sources: the exponential decay of a fixed 
number of ^^Ge atoms and a constant-rate background 
(different for each peak) . Under this assumption the like- 
lihood function |38|| for each peak region is 



N 



where 



/: = e-™]^(6 + ae-^*0, 

i=l 

m = bT + aA/X, 

n 



(6.1) 



fe=i 



e^^*=^). 



fc=i 



Here b is the background rate, A is the decay constant 
of ^^Ge, ti is the time of occurrence of each event with 
t = at the time of extraction, and N is the total number 
of candidate events. The production rate p of ''^Ge is 
related to the parameter a by 



a = ep{l-e-^% 



(6.2) 



where Q — ts — ts is the exposure time (i.e., the time of 
end of exposure Ie minus the time of beginning of expo- 
sure is), and e is the total efficiency for the extraction 
(i.e., the product of extraction and counting efficiencies). 
The total counting live time is given by T and is a sum 
over the n counting intervals, each of which has a start- 
ing time tfjfc and ending time t^k- The parameter A is 



the live time weighted by the exponential decay of ''"'^Ge. 
Its value would be unity if counting began at the end 
of extraction and continued indefinitely. We convert the 
production rate (in ''^Ge atoms produced per day) to the 
solar neutrino capture rate (in SNU) using the conver- 
sion factor 2.977 x 10""* atoms of ^^Ge produced/(SNU 
day ton of gallium) , where the mas s o f gallium exposed 



in each extraction is given in Table HI 



Because of the eccentricity of the Earth's orbit, the 
Earth-Sun distance, and thus the production rate, varies 
shghtly during the year. We correct the production rate 
for this effect by multiplying e by the factor 1 + C where 
C is given by 



C 



with 



2e \ 

— — — ^ ][cosXE+rsmXE 
S[l + J 

-(1 - S){cosXB+rsmXB)], 



r = ti-i/A, 
Xe = ^{tE — tp), 
Xb = t^(is — tp), 



(6.3) 



S = l 



-xe 



Here e is the eccentricity of the orbit (e = 0.0167), ui is 
the angular frequency (w — 27r/365.25 day""'^), and tp is 
the moment of perihelion passage, which has been 2-5 
January for the past number of years. We use tp = 3.5 
days. 

The best estimate of the solar neutrino capture rate 
in each peak region is determined by finding the values 
of a and b which maximize C. In doing so we exclude 
unphysical regions; i.e., we require a > and 6 > 0. The 
uncertainty in the capture rate is found by integrating the 
likelihood function over the background rate to provide 
a likelihood function of signal only, and then locating the 
minimum range in signal which includes 68% of the area 
under that curve. This procedure is done separately for 
the L and K peaks and the results are given in Tables 
IX and ^ We call the set of events in each peak region 
a "data set." 

The overall likelihood function for a single extraction 
is the product of the separate likelihood functions for 
the L- and K-peak regions. The best fit capture rate is 
found by maximizing this function, allowing the indepen- 
dent background rates in the L and K peaks to be free 
variables. The uncertainty in this result is determined by 
finding the values of the capture rate at which the log- 
arithm of the likelihood function decreases by 0.5, again 
choosing the background rates at these two points to be 
those which maximize the likelihood function. The re- 
sults for all extractions that could be analyzed in both 



peaks are given in Table X] 



The capture rate for each extraction of all runs of 
SAGE is plotted in Fig. ||. These results are derived 
from the K peak plus L peak wherever possible, other- 
wise from the K peak alone. For those readers who may 
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TABLE IX. Results of analysis of L-peak events. Nw'^ is a measure of the goodness of fit. See for its calculation and 
use. The probability values are derived from 1000 simulations and have an uncertainty of ~ 1.5%. The total counting live time 
is 10.379 yr. The entries in columns 2-4 include all analysis time cuts. 





Lead 


Live 




Number of 


Number 


Best 


68% conf. 






Exposure 


time 


time 




candidate 


fit to 


fit 


range 




Probability 


date 


(h) 


(days) 


Delta 


events 


'"-Gc 


(SNU) 


(SNU) 




(%) 


Sep. 92 


29.0 


103.8 


0.811 


7 


4.0 


109 


46- 


174 


0.039 


82 


Oct. 92 


27.3 


96.3 


0.839 


10 


0.0 





Q_ 


61 


0.179 


19 


Nov. 92 


30.7 


66.7 


0.688 


12 


0.0 





Q_ 


62 


0.238 


12 


Dec. 92 


26.7 


47.5 


0.835 


10 


7.4 


153 


51— 


208 


0.055 


65 


Jan. 93 


29.9 


23.4 


0.518 


4 


4.0 


135 


35— 


181 


0.092 


67 


June 93 


33.3 


120.7 


0.699 


9 


1.1 


29 


Q_ 


107 


0.490 


2 


Oct. 93-2 


51.9 


71.6 


0.686 


2 


2.0 


193 


19— 


297 


0.097 


55 


Oct. 93-3 


31.6 


102.7 


0.772 


3 


3.0 


287 


88— 


428 


0.078 


68 


July 94 


45.6 


136.0 


0.782 


10 


2.2 


65 




131 


0.026 


95 


Aug. 94 


32.6 


116.4 


0.838 


20 


0.0 





Q_ 


67 


0.056 


73 


Sep. 94-1 


40.5 


120.0 


0.729 


20 


4.7 


171 


54— 


300 


0.087 


32 


Nov. 94 


30.4 


112.3 


0.660 


10 


2.7 


76 


18- 


143 


0.041 


79 


July 95 


35.5 


110.6 


0.776 


16 


1.2 


35 


Q_ 


104 


0.336 


3 


Aug. 95 


35.2 


108.9 


0.698 


16 


3.9 


113 


42- 


200 


0.095 


28 


Sep. 95 


124.3 


80.4 


0.561 


23 


0.2 


8 


Q_ 


179 


0.160 


23 


Oct. 95 


39.3 


120.7 


0.793 


17 


3.2 


169 


33— 


319 


0.041 


78 


Nov. 95 


37.2 


149.9 


0.759 


19 


8.4 


214 


124- 


310 


0.064 


45 


Dec. 95-2 


78.3 


119.9 


0.530 


22 


0.8 


40 


0- 


174 


0.102 


42 


Jan. 96 


33.9 


141.2 


0.767 


21 


0.0 





0- 


61 


0.065 


66 


May 96 


35.2 


117.8 


0.628 


25 


3.5 


104 


23- 


200 


0.038 


82 


Aug. 96 


32.9 


148.7 


0.790 


20 


5.6 


126 


58- 


204 


0.048 


68 


Oct. 96 


33.6 


155.5 


0.785 


11 


0.0 





0- 


48 


0.119 


39 


Nov. 96 


35.0 


162.5 


0.795 


13 


0.2 


5 


0- 


58 


0.042 


85 


Jan. 97 


34.5 


160.0 


0.816 


16 


1.2 


24 


0- 


68 


0.581 


1 


Mar. 97 


36.3 


160.9 


0.814 


10 


2.6 


45 


9- 


89 


0.126 


17 


Apr. 97 


35.0 


167.4 


0.791 


12 


0.0 





0- 


27 


0.108 


45 


June 97 


35.8 


173.4 


0.797 


16 


4.5 


95 


40- 


161 


0.089 


30 


July 97 


37.0 


140.0 


0.752 


13 


0.7 


14 


0- 


61 


0.238 


10 


Sep. 97 


33.6 


166.6 


0.826 


12 


1.1 


24 


0- 


77 


0.059 


64 


Oct. 97 


34.2 


149.5 


0.780 


19 


4.7 


99 


42- 


167 


0.041 


76 


Dec. 97 


34.4 


137.1 


0.726 


15 


3.1 


69 


18- 


131 


0.045 


73 


Combined (31 data sets) 






433 


64.3 


55 


43- 


68 


0.020 


> 99 



be interested in looking for temporal phenomena, the be- 
ginning time tg and ending time for each run can be 
inferred from the mean exposure date tm and total expo- 
sure time Q given in Table [I] by the relationships 



tE = e-tB. 



(6.4) 



B. Global fits 

The combined likelihood function for any set of extrac- 
tions is the product of the overall likelihood functions for 
each extraction. The best fit capture rate for the set of 
extractions is determined by maximizing this function, 
requiring the production rate per unit mass of Ga to 
be the same for each extraction, and allowing the back- 
ground rates in both the L and K peaks to be different 
for each extraction. The uncertainty is found in the same 
way as for the L- and K-peak combination for a single 
extraction. There are a number of other techniques for 
estimating the uncertainties in addition to the two de- 
scribed and used here. When many runs are combined, 
the likelihood as a function of capture rate approaches a 
Gaussian, and the difference between the results of these 
techniques becomes slight. 



The results of global fits to our data are given in 

Sec.lvm. 



VII. SYSTEMATIC UNCERTAINTIES 

There are four basic sources of systematic error in 
SAGE: uncertainty in the chemical extraction efficiency, 
uncertainty in the counting efficiency, uncertainty due to 
nonsolar neutrino production of ^^Ge (such as by cosmic 
rays), and uncertainty due to nonconstant events which 



mimic ' Ge (such as may be made by Rn). Table XII 



summarizes the results of our consideration of all these 
effects and additional information regarding each of these 
items follows. 



A. Chemical extraction efficiency 

The way in which the chemical extrac tion efficiency is 
determined was described in Sec. 
sources of uncertainty. 



IIIB. There are four 
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TABLE X. 
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FIG. 8. Capture rate for each extraction as a function of time. All uncertainties are only statistical. The symbols 1, 2, and 
3 show the combined result for SAGE I, II, and III, respectively. 



TABLE XI. Results of combined analysis of L-peak and 
iiT-peak events for all 31 runs that could be analyzed in both 
peaks. The uncertainty in the probability is ~ 1.5%. Treating 
the L- and K-peak regions as two separate data sets, the total 
counting live time is 21.282 yr. 





Number of 


Number 


Best 


68% conf. 






Exposure 


candidate 


fit to 


fit 


range 




Prob 


date 


events 


^iQe 


(SNU) 


(SNU) 




(%) 


Sep. 92 


13 


6.0 


79 


44- 


123 


0.097 


25 


Oct. 92 


21 


3.3 


32 


4- 


67 


0.105 


26 


Nov. 92 


28 


4.3 


56 


10- 


111 


0.047 


70 


Doc. 92 


28 


16.8 


168 


115- 


229 


0.057 


53 


.Jan. 93 


17 


10.0 


124 


81- 


177 


0.089 


32 


June 93 


18 


3.3 


42 


4- 


92 


0.557 


< 1 


Oct. 93-2 


4 


3.0 


141 


60- 


245 


0.049 


83 


Oct. 93-3 


7 


4.0 


185 


80- 


303 


0.052 


77 


July 94 


22 


3.4 


47 


9- 


94 


0.027 


95 


Aug. 94 


27 


3.9 


46 


15- 


85 


0.075 


52 


Sep. 94-1 


30 


6.5 


112 


50- 


188 


0.082 


39 


Nov. 94 


23 


8.0 


116 


66- 


176 


0.015 


> 99 


July 95 


33 


5.0 


74 


19- 


138 


0.063 


55 


Aug. 95 


24 


7.4 


105 


60- 


161 


0.061 


56 


Sep. 95 


33 


1.2 


28 


0- 


142 


0.058 


73 


Oct. 95 


26 


6.5 


163 


75- 


270 


0.019 


> 99 


Nov. 95 


32 


10.2 


127 


78- 


185 


0.032 


88 


Dec. 95-2 


40 


0.5 


12 


0- 


95 


0.068 


62 


Jan. 96 


35 


3.5 


45 


0- 


101 


0.047 


76 


May 96 


31 


5.3 


78 


31- 


136 


0.039 


90 


Aug. 96 


26 


4.5 


51 


14- 


96 


0.089 


35 


Oct. 96 


21 


5.4 


58 


28- 


95 


0.046 


74 


Nov. 96 


28 


1.9 


21 


0- 


57 


0.103 


37 


Jan. 97 


24 


2.6 


26 


0- 


60 


0.190 


13 


Mar. 97 


24 


6.1 


54 


24- 


90 


0.134 


15 


Apr. 97 


22 


2.7 


27 


3- 


57 


0.037 


86 


June 97 


27 


10.4 


109 


71- 


155 


0.078 


35 


July 97 


22 


0.0 





0- 


24 


0.333 


7 


Sep. 97 


17 


4.3 


49 


22- 


84 


0.043 


80 


Oct. 97 


26 


3.4 


36 


9- 


72 


0.083 


49 


Dec. 97 


24 


6.2 


80 


40- 


128 


0.031 


89 


Combined 


753 


152.1 


64 


56- 


72 


0.033 


93 



1. Mass of Ge carrier 

The extraction efficiency is measured by adding to the 
Ga metal several slugs of Ga-Ge alloy which contain a 
known mass of Ge. This alloy is produced in large batches 
by reduction of Ge by Ga metal from chloride solution, 
and then divided into several hundred small slugs, each 
of which weighs 18-20 g and contains about 40 //g of Gc. 
The equality of Ge content was measured by extracting 
the Ge from a few dozen slugs. The standard deviation 
of these measurements was 2.1%, which we take as the 
uncertainty in the mass of added Ge carrier. 

2. Mass of extracted Ge 

There is also an uncertainty in how much carrier has 
been synthesized into GCH4. This is determined by the 
accuracy to which the GCH4 volume can be determined 
and is estimated to be 2.5%. 



3. Residual Ge carrier 

Since the extraction of carrier Ge is not complete, 
residual Ge carrier from preceding extractions will con- 
tribute to the extraction efficiency measurement. Each 
extraction for solar neutrino data is followed by a second 
extraction to remove this surplus carrier. The amount 
removed during the first two extractions is at least 95%, 
but is uncertain as described above, which leads to an 
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TABLE XII. Summary of the systematic uncertainties. 
The SNU values for extraction and counting efficiency are 
based on a rate of 67.2 SNU. 

Uncertainty 



1. Volume efficiency 



Origin of uncertainty 


in percent 


in SNU 


Extraction efficiency 






Ge carrier mass 


±2.1% 


± 1.4 


Mass of extracted Ge 


±2.5% 


± 1.7 


Residual Ge carrier 


±0.8% 


± 0.5 


Ga mass 


±0.3% 


± 0.2 


Total (extraction) 


±3.4% 


± 2.3 


Counting efficiency 






Volume efficiency 


±1.4% 


± 0.9 


End losses 


±0.5% 


± 0.3 


Monte Carlo interpolation 


±1.0% 


± 0.7 


Shifts of gain 


-3.1% 


±2.1 


Resolution 


±0.5%, -0.7% - 


0.3,±0.5 


Rise time limits 


±1.0% 


± 0.7 


Lead and exposure times 


±0.8% 


± 0.5 


Total (counting) 


±2.3% -3.9% - 


1.5,±2.6 


Nonsolar neutrino production of '^^Ge 




Fast neutrons 




< - 0.02 


232 rjij^ 




< - 0.04 


22«Ra 




< - 0.7 


Cosmic-ray muons 




< - 0.7 


Total (nonsolar) 




< - 1.0 


Background events that mimic 


"Ge 




Internal ^^^Rn 




< - 0.2 


External ^^^Rn 




0.0 


Internal ^^Ge 




< - 0.6 


Total (background events) 




< - 0.6 


Total 




3.0,±3.5 



As described in Sec. [VC, the volume efficiency of 



seven counters of the "LA" type has been directly mea- 
sured with an uncertainty of 0.6%. These seven counters 
were used for 48 of our 88 data sets. The uncertainty in 
El, for counters of this same type used in 36 other data 
sets is estimated from the spread in the measured e„ for 
the measured counters, which is ±2.3% relative uncer- 
tainty. The uncertainty in for counters of the "Ni" 
and "RD" types, which were used for four data sets, is 
taken as ±3%. Averaging over the different types of coun- 
ters used for all extractions, the uncertainty assigned to 
volume efficiency is taken as ±1.4%. 



2. End effects 

The reduced electric field near the ends of the counter 
cathode results in a fraction of events that lie outside 
the ±1 FWHM energy windows. Uncertainties in these 
end effects are due to variations in the physical dimen- 
sions of the counters. Based on measurements of various 
"LA" counters, these dimensional differences lead to an 
uncertainty of ±4.1% in the end effect. This gives ±0.5% 
relative uncertainty in the factor (1 — /d). This should 
be valid for measurements made with the "LA" counters, 
which were used for most of our data, and this value is 
taken for the entire data set. 



3. Monte Carlo interpolation of measured gas parameters 



uncertainty in the amount of remaining carrier. The ex- 
traction efficiency uncertainty due to the uncertainty in 
the residual carrier is ±0.8%. 



4- Mass of Ga 

The total mass of Ga has been weighed periodically 
with a precision of 0.3%. The amount removed during 
each extraction is small (typically 0.1%) and is known 
well (2%). We take the uncertainty in the Ga mass for 
all runs to be ±0.3%. 



B. Counting efficiency 



4.1 



The counter efficiency is calculated using Eq 
There are thus three sources of uncertainty: the volume 
efficiency , the end effects (or equivalently the fraction 
of degraded events), and the gas efficiency. 



The uncertainties in the gas efficiency consist of three 
components: uncertainty in the Monte Carlo calcula- 
tions, uncertainty in the measured gas pressure, and un- 
certainty in the measured percentage of GeH4. The lim- 
ited statistics used in the Monte Garlo calculations to de- 
termine the constants in the gas efficiency formula leads 
to an uncertainty of 1.0% in the determination of the gas 
efficiency. The uncertainty in the gas pressure measure- 
ments is ±5 Torr, which corresponds to an uncertainty in 
the gas efficiency for a typical counter filling (710 Torr at 
24% GeH4) of ±0.2% relative change. The uncertainty 
in the measured percentage of GeH4 is taken to be ±1%, 
which corresponds to an uncertainty in the gas efficiency 
for an average counter filling of ±0.2% relative change. 
Adding these three contributions in quadrature yields a 
relative total uncertainty in the gas efficiency of ±1.0%. 



4. Gain shifts 

If the calibration mean shifts between two calibrations, 
there is an error made in the efficiency estimate. This 
error has been minimized by two features of our stan- 
dard analysis. First, we use a two FWHM wide energy 
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window. Since the peak is relatively Gaussian and the 
window limits are far out on the tail, uncertainties in the 
location of the centroid of the peak do not greatly af- 
fect the efficiency. Second, by sliding the energy window 
between calibrations we hope to minimize any error in 
estimating the centroid due to the observed gain shifts. 
Although the correction for nonlinearity of the counter 
response [Eq. (4.2)] results in an additional uncertainty 
in the gain of 0.7%, the total uncertainty in the gain is 
dominated by the shifts. 

To estimate the error generated by using an incor- 
rect centroid, we computed the area under a Gaussian 
between two integration limits which are shifted by an 
amount 5. We then compared this number to the 0.9815 
number expected from integration limits of ±2 FWHM. 
Using a typical K-peak resolution of 20%"23% we calcu- 
lated the true efficiency for various values of S expressed 
as a fraction of the true mean. 

Typical gain shifts are of the order of a few percent. 
This results in an uncertainty of approximately —3.1% 
in the efficiency. Note that this effect can only decrease 
our efficiency so it is a one sided systematic uncertainty. 



5. Energy resolution 

As a result of the statistics of our calibration spectra, 
the resolution is known to about 2.1%. For the K peak, 
there is an additional uncertainty due to the counter non- 
linearity [Eq. (O)] of ±4.5%. Adding these in quadra- 
ture, the uncertainty in the resolution results in an uncer- 
tainty in the efficiency of about -1-0.5%, —0.7%. Again, 
because the energy window is so wide, the uncertainty 
in the efficiency due to the resolution uncertainty is not 
large. 



6. Rise time limits 



As described in Section |V Q , when the wave form 
method of rise time determination is used, there is an 
uncertainty in the efficiency of ±1% that arises from 
changes in the rise time limits due to counting channel 
and filling variations. For those runs that used the ADP 
method of rise time determination, we can find the un- 
certainty of the ADP cut as follows: Usually a calibra- 
tion has between 1000 and 5000 events in the peak. We 
base our lower ADP threshold on 4% of those or 40-200 
events. This small number of events is subject to sta- 
tistical fluctuations. For the most extreme case, we take 
the square root of 40 (6.4) and notice that the efficiency 
due to the ADP cut could actually be between 94.4% and 
95.6% instead of the 95% we believe it to be. Thus this 
is a ±0.6% uncertainty. Since the vast majority of our 
data are based on wave form analysis, we use ±1% for 
all runs. 



7. Lead and exposure times 

Because extraction usually occurs from several reac- 
tors over the course of 6-10 h, there is an uncertainty 
in the exposure time and in the time from extraction to 
the start of counting (which we call the "lead time" ) of 
roughly 3 h. The lead time is typically 36 h and the expo- 
sure time is typically 34 days. These small uncertainties 
make a small contribution to the uncertainty associated 
with the solar neutrino flux. By Eq. (6.2), the solar neu- 
trino production rate p is proportional to the quantity 



-Ati, 



'(1 



-Ae 



)] ^, where A is the ''^Ge decay con- 



stant, ticad is the lead time, and Q is the exposure time. 
By differentiation one finds that Sp/p due to iiead is about 
±0.8% and due to 9 is about ±0.11%. 



C. Nonsolar neutrino contributions to the '^^Ge 
signal 

In addition to solar neutrinos, ''^Ge can also be pro- 
duced from Ga by the reaction '^^Ga{p, 7T.)^^Ge. The pro- 
tons that initiate this reaction can be secondaries made 
by the (n,p) reaction of fast neutrons or by the {a,p) 
reaction where the a's are from radioactive decay or may 
arise from photonuclear reactions initiated by cosmic-ray 
muons. The yields of these reactions have been mea- 
sured with neutrons from radioactive sources, a's from 
a Van de Graaff generator, and high-energy muons from 
accelerators (see, e.g., Q). Based on these results, great 
care was taken in the design and construction of SAGE 
to minimize these potential background sources. A ma- 
jor advantage of using Ga metal as the solar neutrino 
target (as opposed to an aqueous solution, such as the 
GaCla target of GALLEX) is that the target contains no 
free protons, and thus the production rates of all these 
reactions are low. 

Any one of these processes could produce a background 
effect that must be subtracted from our measured solar 
neutrino signal, but as will be seen below, our best es- 
timates for all of these effects are very small and have 
large errors. Thus, rather than making a background 
subtraction, we include these effects here as systematic 
uncertainties. 

Other Ge isotopes that may be misidentified as ^^Ge 
can be produced in similar reactions: ®^Ge can be made 
by ^^Ga(p, n)^^Ge and the spallation reaction on Ga by 
throughgoing cosmic-ray muons can make ^^Ge and ^^Ge. 
Since the production rate of ^^Ge by the spallation reac- 
tion is comparable to that of ''^Ge and its half- life is long 
(271 days), the ^*Ge decay rate is much less than that of 
^^Ge and can be neglected. The short-lived isotope ^^Ge 
has a greater potential to give events that mimic ^"'^Ge 
and is considered below in Sec. VII D 2 . 
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1. External neutrons 

Fast neutrons mainly arise from the walls of the Ga 
chamber by (a, n) reactions where the a particles are 
from 238^ decay. This background is expected to be small 
due to the low cross sections for (n,p) reactions on Ga 
isotopes and to the low fast neutron background, which 
is because the Ga chamber is lined with low-radioactivity 
concrete and steel. The fast neutron flux in the gallium 
chamber was measured by extracting ^^Ar from a tank 
with 187 kg of dry CaC2 04 and counting in a propor- 
tional counter. Fast neutrons above 3 MeV produce ^^Ar 
through the '*°Ca(n, a)'^^Ar reaction. The flux over 3 
MeV is (4.6 ± 1.6) x lO'^ neutrons/(cm2 day) @. The 
number of ^^Ge atoms produced by this flux in 60 tons 
of Ga metal is < 2.9 x IQ-^'/day |l|-||, which, using 
the conversion factor 56 SNU/(^^Ge atom produced per 
day in 60 tons of Ga), corresponds to < 0.016 SNU. 



2. Internal radioactivity 

The second possible background source is due to a ra- 
dioactivity in the gallium. The only appreciable sources 
of high-energy a's are from decays in the U and Th 
chains, mainly the 8.8-MeV a from ^^^Po at the end of 
the Th chain and the 7.8-MeV a from ^^"^Po near the 
end of the U chain. The concentrations of radioactive 
impurities in the Ga have been measured in two ways: 
by direct counting in a low-background Ge detector by 
the Institute for Nuclear Research (INR) [Q and by glow 
discharge mass spectrometry by both Charles Evans As- 
sociates lH] and Shiva Technologies [|6|. No U, Th, 
or Ra was detected in any of these measurements. Ex- 
pressed in grams of impurity per gram of Ga, the limits 
are U < 2.0 x 10"^° (Evans) and < 1.2 x 10"!° (Shiva); 
Th < 8.0 X 10-1° (INR), < 1.7 X 10^1° (Evans), and 
< 1.2 X 10-1° (Shiva); 226Ra < 1.1 x lO'^^ (INR). We 
take the INR limit for 226Ra and the Shiva limit for Th. 
Using the measured ''^Ge yields |^ in metallic Ga, the 
number of ''^Ge atoms produced per day in 60 tons of 
Ga is < 0.001/day from "^^"^Th and < 0.013/day from 
226Ra, which correspond to < 0.04 SNU and < 0.7 SNU, 
respectively. 



3. Cosmic-ray muons 

The third possible background source is production of 
Ge isotopes by cosmic-ray muons. The global cosmic-ray 
muon flux in the gallium laboratory at BNO has been 
measured Q to be (3.03 ±0.10) x 10"^ muons/(cm2 s). 
This flux can be converted to a ^^Ge production rate in 
60 tons of Ga metal in two ways: (1) using cross sec- 
tions measured at accelerators ||2|,Q and scaling by the 
average muon energy (which is ~ 381 GeV) to the 0.73 
power 1^^, one predicts 0.012/day or (2) using cross 



sections for production of ^^Ge calculated in one 
predicts 0.013/day. Both of these convert to rates of 0.7 
SNU. Since the error of these estimations is about 100% 
p3| , we consider this background as a systematic uncer- 
tainty for the muon background rate. 



D. Background events that mimic '^^Ge 



As is evident from Tables IX and ^ a large fraction 
of the events that we select as^^Ge candidates are not 
''^Ge. In the L [K) peak we select 433 (604) events as 
candidates, but most of these events occur late in the 
counting, and thus the best fit to ^^Ge plus constant 
background assigns only 64.3 (143.7) of them to be ''^Ge. 
These late events that we cannot discriminate from ^^Ge 
are produced by background processes, such as (3 rays 
whose path through the counter either is very short or is 
parallel to the anode wire. As long as these background 
events occur at a constant rate, they only deteriorate 
our signal-to-background ratio, but do not change our 
extracted ^^Ge signal rate. If these background events 
mainly occur early or late in the counting, however, the 
extracted ''^Ge signal rate will be incorrect, too high or 
too low, respectively. Particularly insidious in this regard 
is the ubiquitous naturally occurring isotope ^^^Rn. 



1. Internal radon 

Because it has a short half-life of only 3.8 days, and can 
produce events that mimic ^^Ge, any ^^^Rn that enters 
the counter at the time it is filled will produce events 
early in the counting period that may be falsely inter- 
preted as ^^Ge, and thus give an incorrectly high ^^Ge 
signal rate. To understand this process, let us flrst con- 
sider the principal ^^^Rn decay sequence, which is 

3-82 d I 

3.05 m I 

214Bi_ 



26.8 m 



19.7 m 



Pi 

2i4Po + /32 

L, 210pb_ 



10-'' 
71/ 



"3- 



Events that are falsely identified as ' Ge can be pro- 
duced by one of the a's (this happens rarely as they are 
heavily ionizing), by one of the /3's (this occurs more fre- 
quently), by the recoil nucleus from the a decay if the 
initial nucleus is on the counter wall, or by a low-energy 
X ray emitted by one of the heavy elements in the chain. 
Fortunately, the start of this chain is easily recognized 
as at least one of the first two heavily ionizing a parti- 
cles usually produces a pulse that saturates the energy 
scale. Thus, since this chain takes on average about 1 
hour from the initiating decay of ^^^Rn to reach 22-yr 
^i°Pb, if one makes a time cut of a few hours after each 
saturated event, then most of these false ^^Ge events will 
be removed. 
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TABLE XIII. Rows 1 and 2: probability of a false ^Ge event in the L and K peaks produced by the decay of each nucleus 
in the Rn decay chain per decay of ^^^Rn. Row 3: probability that each nucleus in the chain will decay later than 180 min 
after a ^^^Rn decay. 







Decaying 


nucleus in the ^^^Rn chain 






222Rn 


218pQ 


214p^ 


2143. 


2i*Po 


Prob. false L event 





0.0042 


0.0110 


0.0072 


0.0016 


Prob. false K event 


0.00004 


0.0490 


0.0061 


0.0018 


0.0186 


Survival probability 








0.011 


0.035 


0.035 



We choose to eliminate all events that occur from 15 
min before to 180 min after each detected saturation 
event (energy greater than 16 keV). To determine the 
effect of this time cut, we filled a counter with a typ- 
ical mixture of Xe and GeH4 to which ^^^Rn had been 
added and measured it in system 3 under conditions iden- 
tical to those of solar runs. Based on these measurements 
and Monte Carlo modeling, the spectrum of pulses in the 
counter was determined by each element in the Rn chain. 
The probability of a false ^^Ge event can then be d irectly 



calculated and the results are given in Table XIII 



Folding the probability of a false event with the prob- 
ability of survival after the time cut (Table XIII), we 
obtain the probability of observing a false ^^Ge event af- 
ter the time cut to be 0.00043 in the L peak and 0.00078 
in the K peak. The resolution of this counter was better 
than for the average solar neutrino extraction. If we use 
the average resolution, the number of false '^"'^Ge events 
in a typical solar neutrino run that satisfy all our event 
selection criteria divided by the number of detected sat- 
urated events due to ^^^Rn is calculated to be 0.0006 and 
0.0012 for the L and K peaks, respectively. 

To estimate the number of false ^^Ge events that sur- 
vive the time cut, we next calculate how many saturated 
events are present in our data that can be attributed to 
Rn. This is done by taking the data for each run, making 
the usual time cut after shield openings, and selecting all 
saturated events. The time sequence of these events is 
then fit to a decaying component with the 3.82-day half- 
life of ^^^Rn plus a constant background. For the peri- 
ods of SAGE II and III, this yields 294 saturated events 
initiated by 222Rn, with 192.1 (294.0) events in the 31 
(57) data sets that give the L- {K-) peak results. Since 
SAGE I did not have the capability to detect saturated 
events, we scale the number for the K peak by the num- 
ber of additional extractions (16) to make the iiT-peak 
total 376.5. 

Combining these results, the number of false ^"'^Ge 
events that remain after the time cut is then given by 
0.0006 X 192.1 = 0.11 in the L peak and 0.0012 x 376.5 = 
0.44 in the K peak. Since we have observed a total of 
64.3 (143.7) events in the L [K) peaks, the fraction of 
false events is 0.2% (0.3%), which translates to a false 
signal rate of 0.1 (0.2) SNU. Combining the L- and K- 
peak results gives a total false signal rate of 0.2 SNU. 
Because the uncertainty in this correction is comparable 
with the magnitude of the correction itself, we choose to 



treat this effect as a systematic uncertainty, rather than 
as a background to be subtracted from the signal. 



2. External radon 

^^■^Rn that is external to a proportional counter can 
also produce false ^^Ge events. The radon levels in the 
counting room vary with external conditions, and usu- 
ally fall within the range of (2.0 ± 1.0) pCi/1. To reduce 
the level of Rn in the vicinity of the counters, all passive 
shields are equipped with purge lines from evaporating 
liquid nitrogen and the shields have been made fairly 
hermetic. Each time a counter is calibrated, however, 
the shield must be opened, and some mine air will enter 
the volume around the counters. Under normal circum- 
stances calibrations occur regularly every 2 weeks. Any 
false ^^Ge events that are produced by external Rn will 
thus occur more or less constantly in time, and will be 
treated by the maximum likelihood analysis as a constant 
background. Nevertheless, we minimize the effect of ex- 
ternal Rn by making a time cut on the data for 2.6 h 
after any shield opening. 

A special counter was constructed to give information 
on the false ^^Ge events that are produced by external 
^^^Rn. This consisted of one of our usual counters en- 
closed within a cylindrical quartz capsule (20 mm diam- 
eter). The sealed volume of the capsule was filled with 
air to which ^^^Rn was added to make the total activity 
3.5 nCi. This counter was measured in counting system 6 
which uses the ADP method of rise time determination. 
In the K peak of ^^Ge, after cuts for energy, ADP, and 
Nal, the measured count rate was 0.81 ±0.11 events/min. 

Because the internal volume of the proportional 
counter is shielded by the Fe cathode of 1/3 mm thick- 
ness, the false ^^Ge events are mainly produced by the 
/3 particles from the decay of ^^"^Bi which have sufficient 
energy to penetrate to the active volume of the counter. 
We use the measurements from the internal Rn section, 
our Monte Carlo model for the counter response, make 
some reasonable assumptions regarding the location of 
Rn-daughter products, take into account the reduction 
of Rn in the vicinity of the counters due to the N2 purge, 
and calculate the number of false ^-'^Ge events to be 0.005 
in the sum of the L and K peaks per run. These calcu- 
lations were made without taking into account the effect 
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of the 2.6-h time cut after each shield opening, which 
reduces the number of events even further. Thus the 
number of false ^^Ge events produced by external ^^^Rn 
is negligible. 

This conclusion is verified by analyzing our full L + 
K data set without making the shield opening time cut. 
The result is 67.9 SNU, nearly equal to the result of 67.2 
SNU when the time cut is used. 



3. Internal ^'-^Ge 

Because it can be produced by the same background 
reactions that make ''^Ge, has a short half- life of 39 h, 
and 64% of its decays are by electron capture, ^^Ge can 
produce events that will be misidentified as ''^Ge. We 
can estimate the production rate of ^^Ge from known 
data. The cosmic-ray production rate of ^^Ge can be de- 



termined in the same way as was done in Sec. VII C 3 for 
^^Ge. Using the measured muon flux in the laboratory, 
the cross section for production of ^^Ge of 100 /ib mea- 
sured in GaCla at CERN with 280 GeV muons Q, the 
factor of 2 greater production rate of ^^Ge in Ga metal 
compared to GaCls measured at FNAL with 225 GeV 
muons and scaling as the muon energy to the 0.73 
power, we estimate a production rate of 0.036 atoms of 
^^Ge per day in 60 tons of Ga. The production rate of 
^^Ge by a particles and neutrons is comparable to that 
of ''^Ge, viz., 0.015 atoms/day. We must add to this the 
production rate of ^^Ge by '^B neutrinos, which we esti- 
mate as being comparable to that of ^^Ge, i.e., 5.8 SNU, 
thus making a total estimated production rate of 0.21 
^^Ge/day in 60 tons of Ga. Since the usual exposure in- 
terval is at least 30 days, ^^Ge will be fully saturated, and 
the total number of atoms at the end of exposure will be 
approximately 0.5. When counting starts, on the average 
36 h after extraction, nearly half of these atoms will have 
decayed. Fortunately 86% of the decays of ^^Ge have a 
coincident (3^ or 7 and will be vetoed by the surrounding 
Nal detector with approximately 90% efficiency. Includ- 
ing the 14% of ^^Ge decays that occur by electron capture 
to the ground state, the total efficiency for ^^Ge detection 
will be no more than 25%. Approximately 70% of these 
decays will appear in the L and K peaks, leaving a total 
of 0.045 observed ^^Ge decays per run, or 1 event in every 
44 data sets. Since the 211 events that we have assigned 
to ''^Ge in our 88 data sets correspond to 67 SNU, this 
implies that the false ^^Ge background is approximately 
0.6 SNU. This very small value illustrates the desirabil- 
ity of siting the SAGE detector at great depth and the 
advantage of a Nal veto on all channels during counting. 



VIII. RESULTS 

If we combine SAGE I with SAGE II (minus part 2) 
and SAGE III, the global best fit capture rate for the 88 



separate counting sets is 67.2^7 SNU, where the uncer- 
tainty is statistical only. In the windows that define the L 
and K peaks there are 1037 counts with 211.15 assigned 
to ^^Ge (the total counting live time is 28.7 yr). If we 
were to include the data from SAGE II part 2, the overall 
capture rate would decrease by 7.1 SNU. The systematic 
control of the experiment wa s sus pect during the period 
of the gallium theft (see Sec. [IB), and thus we exclude 
that data interval from our result. 

The total systematic uncertainty is determined by 
ad ding in quadrature all the contributions given in Ta- 
ble KII and is —3.0, -1-3.5 SNU. Our overall result is thus 
67.237^13.0 SNU. If we combine the statistical and sys- 
tematic uncertainties in quadrature, the result is 67.2l7'g 
SNU. 

This section continues with the evidence that we are 
truly counting '^^Ge, considers how well the observed data 
fit the models that are assumed in analysis, and con- 
cludes with consideration of the internal consistency of 
the SAGE results. 



A. Evidence for Ge 

The most direct visual evidence that we are really ob- 
serving ''^Ge is in Fig. ||. The expected location of the 
^^Ge L and K peaks is shown darkened in this figure. 
These peaks are apparent in the upper panel, but miss- 
ing in the lower panel because the ''^Ge has decayed away. 
Events outside the two peak regions occur at about the 
same rate in both panels because they are mainly pro- 
duced by background processes. 

A quantitative indication that ^^Ge is being counted 
can be obtained by allowing the decay constant during 
counting to be a free variable in the maximum likelihood 
fit, along with the combined production rate and all the 
background rates. The best fit half-life to all selected 
events in both L and K peaks is then 10.5^tiJ days, in 
good agreement with the measured value [ pO[ | of 11.43 
days. 



B. Consistency of the data with analysis hypotheses 

1. Energy and rise time window positions 

To test whether or not the energy and rise time win- 
dows are properly set, the windows can be made wider 
and the data reanalyzed. If the rise time window for ac- 
cepted events is increased by 30%, i.e., from 0-10 ns to 
0-13 ns in the L peak and from 0-18.4 ns to 0-24.0 ns in 
the K peak, then the overall result of all runs of SAGE II 
and III that were counted in system 3 is 68.3 SNU. This 
change is entirely consistent with the ^ 3% increase in 
counting efficiency due to the increased size of the rise 
time acceptance window. Similarly, if the energy win- 
dow in both L and K peaks is opened from the usual 2 
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FWHM to 3 FWHM, then the overall result of all runs of 
SAGE II and III becomes 69.1 SNU. This increase from 
the value of 67.2 SNU in the 2 FWHM energy window is 
because some of the ''^Ge decays occur at the ends of the 
counter and their detected energy is reduced from the full 
peak value. This results in an increase in the counting 
efficiency in the wider energy window of 2%-3%. If this 
efhciency increase is included in the analysis, then the 
results in the two energy windows agree to better than 
1%. 



2. Time sequence 

A major analysis hypothesis is that the time sequence 
of observed events for each run consists of the superposi- 
tion of events from the decay of a fixed number of ^^Ge 
atoms plus background events which occur at a constant 
rate. The quantity Nv?' and the goodness of fit proba- 
bility inferred from it provide a quantitative measure of 
how well the data fit this hypothesis (see for the def- 
inition and interpretation of Nw'^). These numbers are 
evaluated for each data set and are given in Tables 1^, 
and XI . There are occasional runs with rather low prob- 



ability of occurrence, but no more of these are observed 
than are expected due to normal statistical variation. 

This method can also be used to determine the good- 
ness of fit of the time sequence for any combination of 
runs. These numbers are given in the various tables; for 
the combined time sequence of all L plus K events from 
all runs, this test yields Nw'^ — 0.074, with a goodness- 
of-fit probability of (58 ±5)%. A visual indication of the 
quality of this fit is provided in Fig. |^ which shows the 
count rate for all events in the L and K peaks vs time 
after extraction. As is apparent, the observed rate fits 
the hypothesis quite well. 



3. Production rate sequence 

Another analysis hypothesis is that the rate of ^^Ge 
production is constant in time. By examination of Fig. H, 
it is apparent that, within the large statistical uncer- 
tainty for each run, there are no substantial long-term 
deviations from constancy. 

To quantitatively test whether or not it is reasonable 
to assume that the production rate is constant, we can 
consider the three segm ents of SAGE data, whose results 
are given in Table XIV. A test of the consistency of any 
data segment with the overall result of 67 SNU can be 
made by Monte Carlo simulation. For the purposes of 
illustration, we choose the most deviant segment, SAGE 
III, whose overall result is 56 SNU. We then simulate all 
39 data sets of SAGE III assuming that the true pro- 
duction rate is 67 SNU. To ensure that these simulations 
parallel the real data as closely as possible, all parameters 
of the simulation, such as background rates, efficiencies, 
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FIG. 9. Gount rate for all runs in L and K peaks. The 
solid line is a fit to the data points with the 11.4-day half-life 
of ^'^Ge plus a constant background. The vertical error bar on 
each point is proportional to the square root of the number of 
counts and is shown only to give the scale of the error. The 
horizontal error bar is ±5 days, equal to the 10-day bin size. 



exposure times, and counting times, are chosen to be the 
same as for the real data. From the sequence of simulated 
event times, the combined production rate is calculated 
in exactly the same manner as for the real data. This 
process is repeated 10000 times and a histogram of the 
combined rate is produced. From the position of the ob- 
served rate for the real data in this histogram, we can 
calculate the probability that the real data are produced 
by the assumed initial production rate. As shown in Fig. 
0, we find that (11 ± 0.3)% of the 10000 simulations of 
SAGE III have a value that is lower than the observed 
value of 56 SNU. Since this probability is one tailed (max- 
imum of 50%), this is the most aberrant of the three 
sections of SAGE data, and no systematic uncertainties 
were included in the simulations, a value of 11% is not 
extremely unusual, and there is thus no statistically sig- 
nificant evidence for production rate variation. The same 
analysis applied to SAGE I and SAGE II yields proba- 
bilities of 35% and 38%, respectively, highly consistent 
with the assumption of constant production rate. 

Another way to consider this question is to use the 
cumulative distribution function of the production rate 
C{p), defined as the fraction of data sets whose produc- 
tion rate is less than p. Figure |ll| shows this distribution 
for all data sets and the expected distribution from sim- 
ulation, assuming a constant production rate of 67 SNU. 
The two spectra parallel each other closely and can be 
compared by calculating the Nw'^ test statistic |^ . This 
gives Nw'^ = 0.343 whose probability is 10%. 
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TABLE XIV. Results of combined analysis of various segments of SAGE data. The time intervals for each segment are 



defined in Table |n The uncertainty in the probability is ~ 4%. 
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FIG. 10. Distribution of capture rate from 10000 simula- 
tions of SAGE III assuming true production rate of 67.2 SNU. 
The probability of a rate less than or equal to the observed 
rate of 56 SNU is 11% and is shown shaded. 
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(jagged curve) and the expected distribution derived by 1000 
Monte Carlo simulations of each set (smooth curve). The 
capture rate in the simulations was assumed to be 67.2 SNU. 



Although these statistical tests are consistent with a 
constant production rate, they can never exclude the 
possibility of a cyclic time variation whose magnitude 
is comparable wit h the statistical uncertainty. We thus 
give in Table KV the capture rate result for several of 
the possible temporal combinations of SAGE data. Each 
of these data divisions fits well to the constant rate of 
67 SNU, as is verified by x^/degree of freedom = 8.2/7 
(yearly), 14.6/11 (monthly), 4.9/5 (January + February 
bimonthly), and 3.9/5 (February + March bimonthly), 
which have probabilities of 32%, 20%, 43%, and 56%, re- 
spectively. We remind those readers who are interested 
in short-term periodicity that the known variation due 
to the change in Earth-Sun distance ha s be en removed 
from our reported capture rate [see Eq. (6.3)]. 



is 6 SNU above. The statistical Itr error of these results, 
however, extends upward to 68 SNU in the L peak and 
downward to 64 SNU in the K peak. Both L- and K- 
peak results thus overlap the overall value, and there is 
no evidence for inconsistency between the results in the 
L and K peaks. 
As noted in Sec. 



[II B 



so as to remove most of the 
residual Ge carrier from the Ga metal, it is customary to 
make a second extraction 2 or 3 days after each solar neu- 
trino extraction. Although these second extractions are 
usually counted, until recently they were often measured 
in counters which did not have the lowest background 
rates, and were rarely counted in electronic system 3 with 
the wave form recorder. Further, these runs were seldom 
counted for a long time. As a consequence, it was not 
possible for us to give a result for the production rate 
from these second extractions. This situation changed 
at the beginning of 1996, however, because SAGE then 
switched to a 6-week extraction schedule, which freed 
some better low background counters and made it possi- 
ble to measure these samples from second extractions in 
below the overaUvaKf 67 SNu'and the X-peak resuft ^y^^em 3. Ten such extractions have been measured since 



C. Internal consistency of SAGE results 

The combined results for all runs in the L and K peaks 
are given in Table [XIV|. The L-pcak result is 12 SNU 
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TABLE XV. Capture rate results for yearly, monthly, and 
bimonthly combinations of SAGE data. Runs are assigned to 
each time period by their mean exposure time. 
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1996. Taking into account the delay between the first 
and second extractions and the extraction and counting 
efficiencies, in these ten extractions we expect to detect 
three ''^Ge atoms that are leftover from the first extrac- 
tion, and seven ^^Ge atoms that are produced by solar 
neutrinos during the interval between extractions. The 
total number of ^^Ge atoms detected in these ten extrac- 
tions was 1.1 with a 68% confidence range from 0.0 to 
8.7. The number observed is statistically consistent with 
the number expected, thus confirming our extraction ef- 
ficiency. Further, it establishes that the ^^Ge we detect 
is not an artifact of the extraction process and that our 
counting and data analysis do not find a significant quan- 
tity of ^^Ge if it is not present. 



IX. SUMMARY AND CONCLUSIONS 

We have presented the methods and procedures of the 
SAGE experiment: the extraction of Ge from Ga, the 
subsequent Ge purification, the counting of ''^Ge, the 
identification of candidate ^^Ge events, and the analysis 
of the counting data to obtain the solar neutrino pro- 
duction rate. Eight years of measurement of the solar 
neutrino flux give the capture rate result 67.2t7 Q SNU, 



where the uncertainty is statistical only. Analysis of all 
known systematic effects indicates that the total system- 
atic uncertainty is g SNU, considerably smaller than 
the statistical uncertainty. Finally, we have examined 
the counting data and shown that there is good evidence 
that ^^Ge is being counted, that the counting data fit 
the analysis hypotheses, and that the counting data are 
self-consistent. 

The SAGE result of 67.2 SNU represents from 52% 
iQ to 53% [H of SSM predictions. Given the extensive 
systematic checks and auxiliary measurements that have 
been performed, especially the ^^Cr neutrino source ex- 
periment | |20| , pT[ |, this 7a reduction in the solar neutrino 
flux compared to SSM predictions is very strong evidence 
that the solar neutrino spectrum below 2 MeV is signifl- 
cantly depleted, as was previously shown for the flux 
by the CI and Kamiokande experiments. If we take into 
account the results of all experiments, astrophysical so- 
lutions to the solar neutrino deflcit can now nearly be 
excluded |^-^. This conclusion is indeed implied by 
the SAGE result itself, as it lies 2.5a below the capture 
rate prediction of 88.1^2 4 SNU obtained by artiflcially 
setting the rate of the '^He(a, 7)^Be reaction to zero and 
1.5a below the astrophysical minimum capture rate of 
79.5^2:0 SNU 0. The solar neutrino problem is now 
a model-independent discrepancy |l6|jr^ ] that does not 
depend on the details of solar models or their inputs. 

More credible explanations for the solar neu- 
trino deficit involve either matter-enhanced Mikheyev- 
Smirnov-Wolfenstein (MSW) neutrino oscillations, in 
which the solar oscillates into other flavor neutrinos 
or a sterile neutrino [p^ , p4| -[56| , or vacuum oscillations 
p7| , p^ 58 1 . For both of these possibilities, the allowed re- 
gions of Am^ — sin^ 29 parameter space determined from 
solar neutrino experiments for two-flavor oscillations into 
active neutrino species are shown in Fig. |l2|. The flt qual- 
ity is about the same in both regions. There is also a 
flt with similar quality for MSW oscillations into sterile 
neutrinos, whose allowed region approximately coincides 
with the region shown for MSW oscillations with active 
neutrinos. 

There are now very strong indications that the solar 
neutrino deficit has a particle physics explanation and is a 
consequence of neutrino mass. To fully unravel the solar 
neutrino story, however, will require more experiments, 
especially those with sensitivity to low-energy neutrinos 
or to neutrino flavor. SAGE continues to perform reg- 
ular solar neutrino extractions every 6 weeks with ^ 50 
tons of Ga and will continue to reduce its statistical and 
systematic uncertainties, thus further limiting possible 
solutions to the solar neutrino problem. 
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APPENDIX A: OTHER COUNTING SYSTEMS 

The counting systems have been designated by the 
numbers 1-6. The initial developmental work on sys- 
tem 1 [ ^9p0[ , which used the amplitude of the differen- 
tiated pulse (ADP) method |6^ to separate ^^Ge events 
from background, was done in Russia during the early 
1980s. Based on this work, system 2 was developed at 
BNO during the years 1985-1988. System 2 was com- 
pleted in 1989 and counted all but two first extractions 
through May 1992 (SAGE I). Counting system 5, which 
used the ADP method of rise time measurement, was 
used to count the other two first extractions during 1990 
and 1991. During the summer of 1992, system 3, which 
has the capability to record the counter wave form, was 
brought on line; since that time, it has been used to count 
almost all first extractions. After the implementation of 
system 3 as the primary counting system, extensive up- 
grades to reduce backgrounds were performed on system 
2 to enable SAGE to have low-noise counting capability 
in more than eight channels. The upgraded system is 
referred to as system 6. It has counted seven first ex- 
tractions during SAGE II and III, mostly from low-mass 
samples of Ga, and has been used mainly for develop- 
mental work, such as testing proportional counters and 
counting cleanup extractions of gallium. 



a. Counting system 2 

System 2 was a seven-channel system where each PC 
was counted in an independent passive shield; five of 
those channels had active shielding with Nal crystals. 
The passive shield consisted of an internal wall of tung- 
sten (40-80 mm thick) or copper (20-30 mm thick) sur- 
rounded by lead (150 mm thick). The Nal events were 
recorded in coincidence mode with events from the pro- 
portional counters. Several of the per form ance charac- 
teristics of system 2 are given in Table VII, 



h. Counting system 6 

Modifications to system 2 began during 1992 when sys- 
tem 3 became the primary acquisition system; the im- 
provements were so extensive that it was redesignated 
as system 6. The counting system has seven channels 
of acquisition with independent passive shields for each 
proportional counter. Six channels have an active shield, 
which operates in coincidence mode with events in the 
proportional counter. A modified ADP method with 
the application of several differentiation time constants 
is used for rejection of point ionization events from back- 
grounds. System 6 became fully operational in early 
1993. To give this system wave form recording capability 
a digitizing oscilloscope was added, but this improvement 
has never been fully implemented. 
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FIG. 13. Determination of the time and dc offset of a can- 
didate event. 



APPENDIX B: FOURIER TRANSFORM OF THE 
WAVE FORM - R 

In contrast to the Tat method, the pulse offset is de- 
termined independently from the wave form. One uses 
the intersection point of two lines, the zero-slope line of 
the offset and the initial slope of the pulse, to obtain the 
onset position in time and voltage. The initial slope is 
defined as a certain number of points before and after the 
point at 20% of the maximum pulse height. The exact 
number of points to fit is determined individually for each 
pulse since the number of points available will depend on 
the pulse height. This region is chosen so that the points 
are sufficiently linear. Figure ^ illustrates graphically 
how the onset point is determined. 

Two data runs with counters filled with ^^Ge were used 
to check the energy offset from this wave form analysis 
determination. The runs were separated by three years 
and used different digitizer settings. The Gaussian cen- 
troid of each L peak and K peak was calculated, with 
each peak containing a few thousand counts. The extrap- 
olated intercepts in energy are 0.005(6) keV and 0.022(8) 
keV using L- and iiT-peak energies of 1.17 keV and 10.37 
keV. Given the energy resolution of our counters, the en- 
ergy offset is effectively zero. 

The algorithm for determining the pulse onset was 
checked using computer-simulated pulses, both with and 
without Gaussian noise. It correctly identifies the time 
offset to within 1 ns and the dc offset to within one chan- 
nel. Those limits are, of course, dependent on the noise 
levels, but the levels used were approximately the same 
as for typical data. Thus, if each pulse is properly nor- 
malized to both zero time and dc offset, there is no need 
to apply an energy offset correction. 

This technique uses the zero- and lowest-frequency val- 
ues from a FFT to obtain measures of the energy and 
rise time of a pulse. The determination of the energy is 
straightforward from the definition of the Fourier trans- 
form, 



At w = 0, F(uj) equals the area under the curve /(i), 
which in this case is the digitized wave form of the event 
convoluted with a Hanning windowing function. We se- 
lect an integration time of 800 ns, which is the maximum 
time allowable given the variation in the time of pulse 
onset. In effect, this technique is equivalent to summing 
channels used with T-^ and is analogous to an ADC that 
integrates for 800 ns. 

In a Fourier analysis, the rise time behavior of a typi- 
cal pulse clearly will be a very- low-frequency component. 
Studies with actual ^^Ge pulses and computer-simulated 
pulses generated with Eqs. (5.1) show that one can accu- 
rately identify several distinct features of the wave form. 
As expected, the dominant components are the lowest 
frequencies along with the random noise that spans all 
frequencies. One can identify structure as well; most 
of it originates from the intrinsic properties of the os- 
cilloscope, such as dithering and the finite digitization 
size. One of the advantages of a Fourier analysis is that 
such structure appears at high frequencies and is well 
separated from the rise time information. The lowest, 
nonzero, real component F(\) scales similarly to an ADP 
value but is independent of electronic offsets and high- 
frequency noise contributions to the pulse. Dividing it 
by the energy i^(0) of the pulse produces a parameter R 
that is proportional to the inverse rise time. Thus, one 
can perform a complementary analysis of the data that 
is analogous to the ADP method but is based solely on 
the digitized pulse and is independent of any underlying 
assumptions of its functional form. 



APPENDIX C: RST METHOD 

In the standard analysis of our data we use the T/v 

This func- 



method and fit the observed pulse to Eq. 5.1 



tion gives the correct description of the shape of the volt- 
age pulse as recorded by the digital oscilloscope when the 
ionization produced in the proportional counter consists 
of a set of point ionizations evenly distributed along a 
straight track. Since ^^Ge events are usually a single 
cluster of ionization, this method works satisfactorily to 
select ''^Ge candidate events. It is, however, restricted 
to the particular form of ionization that is assumed, and 
gives a poor fit to other types of charge deposit in the 
counter, such as the combination of a point event from 
''^Ge if-electron capture followed by capture of the 9.3- 
keV X ray at some other location in the counter. To give 
us the capability to investigate all possible events that 
may occur in the counter, we have also developed a more 
general method which can analyze an event produced by 
ionization with an arbitrary distribution of charge. We 
call this the "restored pulse method" or "RST method" 
for short. 
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We begin with the measured voltage pulse V{t) as 
recorded by the digitizer. For an ideal point charge that 
arrives at the counter anode wire, V{t) has the Wilkin- 
son form V{t) = W{t) = Voln(l + t/to), provided the 
counter is ideal and the pulse processing electronics has 
infinite bandwidth. For a real event from the counter, 
with unknown charge distribution, V{t) can in general 
be expressed as the convolution of the Wilkinson func- 
tion with a charge collection function G{t): 



V{t) ^W{t)(g)G{t). 



(CI) 



The function G{t) contains within it the desired infor- 
mation about the arrival of charge at the counter anode, 
coupled with any deviations of the counter or electronics 
from ideal response. Equation (CI) can be considered as 
the definition of G{t). 

To get the desired function G{t), one must deconvo- 
lute Eq. (CI). To perform this deconvolution, we have 
found it mathematically convenient to use the current 
pulse I{t), which is obtained by numerical differentiation 

ofv{ty. 
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where W it) is normalized over the observed time of pulse 

measurement, Tobs, such that J^"^' W it)dt — 1. 

To deconvolute, we Fourier transform to the frequency 
domain and then use the theorem that convolution in 
the time domain becomes multiplication in the frequency 
domain |6^. This simply gives /(/) = W (/)G(/), which 
can be solved for G(/). We then Fourier transform G(/) 
back to the time domain to get the desired function G(<). 

The energy of the event is given by Jf^^"^" Git)dt. The 
duration of the collection of ionization is given by the 
width of Git), which can be used as a measure of the rise 
time. 

An example of this procedure as applied to a typical 
^^Ge K-peak event is given in Fig. This pulse has 
Tjyf = 3.9 ns. The recorded voltage pulse after inversion 
and smoothing is given by ViT) in the lower panel. The 
current pulse, obtained by numerical differentiation of 
the voltage pulse, is given by lit) in the upper panel. The 
deduced function G(t) is also shown in the upper panel. 
It has a FWHM of about 15 ns, found to be typical for 
true ^^Ge K-peak events. The integrated current pulse, 
which records the pulse energy, is given by / Git)dt in 
the lower panel. 

This method has the advantage that it can reveal the 
basic nature of the ionization in the counter for an arbi- 
trary pulse. It is also capable of determining the pulse 
energy over a wider range than the Tjv method. A prob- 
lem that has been found with this method in practice, 
however, is that when ^^Ge data are analyzed one ob- 
tains multiple collection functions [i.e.. Git) has several 
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FIG. 14. Analysis of typical '^^Ge pulse by the RST 
method. See text for explanation. 



distinct peaks separated in time] more often than is ex- 
pected from the known physical processes that take place 
in the counter. These multiple peaks are due to noise on 
the pulse and cutoff of the system frequency response at 
about 100 MHz. Attempts have been made to remove 
these extraneous peaks by filtering and smoothing the 
original pulse, but they have not been fully successful. 
Evidently we need faster electronics and a reduction in 
the noise level to be able to fully exploit this pulse shape 
analysis technique. As a result, we have only been able 
to use this method to select events on the basis of energy. 
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